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THREE-DIMENSIONAL TRANSONIC FLOW THEORY 
APPLIED TO SLENDER WINGS AND BODIES 

By Max. A. Heaslet and John R. Spreiter 1 


SUMMARY 


The present paper re-examines the derivation of the integral 
equations for transonic flow around slender wings and bodies of 
revolution , giving special attention to conditions resulting from 
the presence of shock waves and to the reduction of the relations to 
the special forms necessary for the discussion of sonic flow, that 
is, flow at free-stream Mach number 1. In the vicinity of the 
body, the disturbance field is then shown to consist of a two- 
dimensional disturbance field extending laterally and a longi- 
tudinal field that depends on the streamwise growth of cross- 
section area. This result extends Oswatitsch- s equivalence rule 
to lifting cases, provided the angle of attack is small relative to 
the thickness ratio. The correctness of the analysis is checked 
by examination of Yoshihara's numerical solution for sonic 
flow around a slender , circular cone-cylinder and this solution 
is checked , in turn, by comparison with experimental results of 
Solomon. An example is presented in which the general result 
is applied to calculate pressure and integrated forces on a family 
of slender, elliptic cone-cylinders. An expression is derived 
which permits the ready calculation of the difference in drag of 
two slender bodies having the same longitudinal distribution of 
cross-section area. Classes of wings and- bodies are described, 
for which the difference in drag is zero and the Whitcomb area 
ride applies. Experimental data for such a family of wings of 
rectangular plan form are examined and it is shown that theory 
and experiment are in good accord. 

INTRODUCTION 

The equations governing transonic flows are known and 
well established by favorable comparisons with experiment 
(see ref. 1 for a resume). The difficulties arising as a result 
of the nonlinearity and mixed character of the differential 
equation for the potential, however, have prevented the 
rapid advancement of the analysis such as has occurred in 
recent years with both subsonic and supersonic theory. 
This is particularly true for three-dimensional transonic 
Hows and, as a result, perhaps greater than usual effort has 
gone into the determination and utilization of relations 
between solutions. The first of these to be advanced was 
the transonic similarity rule which pertains to the pressures 
and forces on affinely related wings (refs. 2, 3, and 4) and 
bodies of revolution (ref. 5). A second relation is the area 
rule established empirically by Whitcomb (ref. 6) which 


states that “near the speed of sound, the zero-lift drag rise 
of thin low-aspect-ratio wing-body combinations is primarily 
dependent on the axial distribution of cross-sectional area 
normal to the air stream.” A third relation is the equiva- 
lence rule of Oswatitsch (refs. 7 and 8) which may be stated 
as follows: “The solution for transonic flow around a thin, 
nonlifting, low-aspect-ratio wing can be obtained from that 
for a nonlifting body of revolution having the same longi- 
tudinal distribution of cross-sectional area by superposing 
the difference between the two-dimensional harmonic cross- 
flow solutions for the two bodies.” The area rule and the 
equivalence rule are, obviously, closely related. Further 
effort needs to be expended, however, in establishing the 
generality and range of validity of these relations and in 
exploiting the results in specific applications. The present 
paper is concerned with this task. 

Tiie problem will be approached through application of 
the classical method of singularities. This is one of the 
oldest and most fruitful methods for solving partial differ- 
ential equations and has reached a high state of development 
in linearized compressible-flow theory. There is also a con- 
siderable body of literature in which the method is applied 
to nonlinear compressible-flow problems by considering the 
solution of the linearized equations to be a first approxi- 
mation, and iterating to obtain second and higher order 
approximations. The results so calculated are good approxi- 
mations to pure subsonic flows or to pure supersonic flows, 
but it is now generally agreed that the series representation 
of the solution does not converge in the transonic range. 
Approximate calculations by Oswatitsch (refs. 9 and 10) 
have indicated the possibility, however, that the method of 
singularities might be applied successfully in the transonic 
range if the idea is relinquished that the linear solution is 
necessarily the first approximation in the transonic range. 
This idea has been pursued further in references 11 through 
15 in which a number of improvements are introduced and 
the results of numerous specific calculations are shown. 
Although the basic equations are derived for three-dimen- 
sional flow in the latter references, all applications are to 
two-dimensional flows. The values of the free-stream 
Mach number, moreover, are restricted to values no greater 
than unity. 


1 Supersedes NACA TN 3717 by Max. A. Heaslet and John R. Spreiter. 
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The same general approach has been applied to three- 
dimensional transonic flow around slender wings and bodies 
by Oswatitsch and Keune (ref. 8) and by Harder and 
Klunker (ref. 16). In these applications, the principal aim 
is not to determine actual solutions but to derive relations 
between solutions for various bodies having the same longi- 
tudinal distribution of cross-section area. These two anal- 
yses are not entirely satisfactory in a number of particulars, 
not the least of which is the omission of all considerations of 
shock waves in the body of the analysis. A more important, 
although perhaps more subtle, point concerns the treatment 
of the cumulative effect of the nonlinear term of the transonic 
equation when the free-stream Mach number M m is unity. 
Harder and Klunker argue that the effect on the induced 
flow field is small because the term itself is everywhere small. 
Actually, however, the cumulative effect of this term leads 
to infinite contributions at M co = 1. Oswatitsch and Keune 
consider the cumulative effect but circumvent the difficulty 
by introducing rather arbitrarily selected values for Mach 
number so chosen that the value of unity is never inserted 
into the vital integrals. It is the initial concern of the 
present analysis, therefore, to re-examine the derivation of 
the integral equations for transonic flow around slender 
wings and bodies of revolution, giving special attention to 
conditions resulting from the presence of shock waves and 
to the reduction of the relations to the special forms neces- 
sary for sonic flow. In contrast to references 8 through 15, 
which are concerned exclusively with cases in which the 
free-stream Mach number is no greater than unity, equations 
are also derived herein for the case where the free stream is 
supersonic. These equations are likewise reduced to the 
special form associated with sonic flow and the results are 
shown to be identical to those which arise from a considera- 
tion of flows with a subsonic free-stream velocity. 

Following the establishment of the basic integral relations 
for transonic flow, special attention is directed toward the 
case where the free-stream Mach number is unity. Here, 
the integral relations are simpler in character, although still 
nonlinear. Application of a convergent interation process 
leads to the conclusion that the solution for the potential 
has a particularly simple form in the vicinity of the body; 
in common with linearized slender-body theory, the dis- 
turbance field consists, to a given order of error, of a two- 
dimensional disturbance field extending laterally and a 
longitudinal field that depends on the streamwise growth 
of cross-sectional area. This result extends Oswatitsch ’s 
equivalence rule to lifting cases, provided the angle of attack 
is small relative to the thickness ratio. The correctness of 
the analysis is checked by examination of Yoshihara’s numer- 
ical solution for sonic flow around a slender, circular, cone- 
cylinder given in reference 17, and this solution is checked, 
in turn, by comparison with experimental results of Solomon 
given in reference 18. The results yield a simple means of 
determining the pressure distribution on an entire family of 
slender wings and bodies having the same longitudinal dis- 
tribution of cross-sectional area when the pressure distribu- 
tion is known for any member of the family. Starting with 
the known solution for the circular cone-cylinder, an example 
is presented in which the general result is applied to a family 


of slender elliptic cone-cylinders. This example, which was 
discussed briefly in reference 1, is examined in detail. It is 
shown that the lift and the load distribution are the same as 
given by linear theory, confirming the ideas advanced in 
refe rence 19. Contrary to Whitcomb’s area rule, however, 
the drag depends significantly on the cross-section shape. 
Both the drag and lift of the thin elliptic cone-cylinder are 
shown to be in accord with the transonic similarity rules. 
A momentum analysis of the sonic drag of slender bodies in 
general is then undertaken and an expression is derived which 
permits the ready calculation of the difference in drag of 
two slender bodies having the same longitudinal distribution 
of cross-section area. This result confirms the drag variation 
calculated for the elliptic cone-cylinders by integration of 
the surface pressures. Several large and significant classes 
of wings and bodies are described for which the difference 
in drag is zero and the Whitcomb area rule applies without 
modification. One of these is a family of affinely related 
wings. Experimental data from reference 20 for such a 
family of wings of rectangular plan form are examined and 
it is shown that theory and experiment are in good accord, 
provided the product of aspect ratio and cube root of the 
thickness ratio is, in this instance, less than about unity. 

The final section of the report has been written so as to be 
as self-contained as possible and readers concerned solely 
with applications of the theory may find this section sufficient 
for their purposes. The initial sections of the report have 
been written for readers concerned with a more complete 
understanding of the derivation and limitations of the 
general theory together with the evaluation of the order 
of error incurred in the slender-body approximation. 

FUNDAMENTAL EQUATIONS AND BOUNDARY CONDITIONS 

The basic equations necessary for the discussion of in viscid 
transonic flow consist of a set of partial differential equations 
relating the velocity components and their gradients at every 
point, together with an auxiliary relation for the velocity 
jump through a shock wave. For thin wings and slender 
bodies inclined at zero or small angles of attack, the differ- 
ential equations can be simplified if the shock waves are 
assumed sufficiently weak that the flow is irrotational and 
isen tropic, thereby permitting the introduction of a velocity 
potential 4>. The further assumption of small disturbances 
leads to the use of a perturbation velocity potential (p which 
in Cartesian coordinates satisfies the following nonlinear 
partial differential equation 

( 1 — MJ)<P xx + <Pyy + *P 2Z = MJ <P X <P XX ( 1 ) 

C oo 

where U m and M m refer to the velocity and Mach number 
of the undisturbed flow, y is the ratio of specific heats 
(7=1.4 for air), and x , y , and z are Cartesian coordinates. 
The perturbation velocity vector is given by the gradient of 
cp and has components, v , and w along the three axes. 

Knowledge of methods for obtaining solutions of equation 
(1) is meager not only because it is nonlinear, but because it 
changes type (elliptic, parabolic, hyperbolic). This change 
of type is an essential feature of transonic flow, since subsonic 
flows are represented by elliptic equations and supersonic flows 
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by hyperbolic equations. If both types of flow occur in a single 
flow field, it is apparent that the differential equation must 
change type. In the present case, the type of the equation is 
recognized by the sign of the total coefficient of as follows 

C >0 elliptic (subsonic) 

(1 — MJ) —MJ( 7+1) jj~ j =° parabolic (sonic) (2) 

hyperbolic (supersonic) 

In most of the investigations of two-dimensional transonic 
flows (<p yv = 0), the differential equation is transformed into a 
linear equation of mixed type (Tricomi equation) by the in- 
troduction of the hodograph variables. At the present time, 
however, no transformation is known that achieves a corre- 
sponding linearization of the three-dimensional equation, and 
the investigation of other methods of solution thus becomes 
relatively more important. 

Equation (1) is, of course, valid only in regions where the 
necessary derivatives exist and are continuous. Since these 
conditions do not hold where shock waves occur, and since 
shock waves are a prominent feature of most transonic flows, 
an additional equation is needed for the transition through 
the shock. The fundamental properties of a shock surface 
require that the normal component of velocity be discon- 
tinuous and the tangential component, and therefore <p, be 
continuous. The necessary relation follows from the classical 
expression for the shock polar, which in the small disturbance 
transonic theory is approximated by 


(i -M ro 2 )K-^ 2 +K-^ 2 +K-^) 2 




7+1 <£++<£. 


u a 


PK -*0 2 ® 


where the subscripts a and b refer to conditions ahead of the 
behind the shock. 

Equations (1) and (3) are usually developed for the case 
where the coordinate system is placed so that the x axis is 
parallel to the undisturbed stream at infinity, but they also 
apply to tiie case where the coordinate system is rotated 
slightly. In the present analysis, it is convenient to aline 
the x axis with the longitudinal axis of the wing or body as 
shown in figure 1. Such a system is usually referred to as 
the body axes. With these coordinates, the relation be- 




tween the total velocity potential $(x,y,z) and the perturba- 
tion velocity potential (p(x,y,z ) is approximated by 

$(x, y, z ,) =U a>(x+ctz)+<p(x, y, z) (4) 

where a is the angle of attack. 

The expression for the pressure coefficient C v is not 
invariant with respect to small rotations of the coordinate 
system. In body axes, the proper expression is 

C p =—jJ- (Vx+ a< Pz)~-jj-2 + (5) 

The boundary conditions require that the gradient of the 
total velocity potential evaluated infinitely far from the air- 
craft be consistent with the uniform free-stream conditions 
there and, when evaluated normal to and on the surface of 
the airplane itself, be zero. The condition at infinity yields 
$ ( oo )= U m (x +<*s) or that 

<p( °o ) =0 (6) 

An exception to this statement occurs in the vicinity of the 
wake at great distances behind the wing, but no complica- 
tion ensues due to the relative smallness of this region. 
The condition at the airplane surface results in the relation 

^^UJni+cxn^+ri! ^ ~+n 2 ~+«- 3 |f=0 (7) 

where n u n 2 , and n 3 are the direction cosines of a normal to 
the airplane surface with respect to the x, y, and z axes, 
respectively. This relation is too general for the present 
needs, however, because it applies to all shapes, whereas the 
analysis is to be a small disturbance theory that applies 
only to slender bodies and thin wings. For such configura- 
tions, rii is small nearly everywhere on the surface and will 
be neglected in comparison with either unity or {n 2 2J rn z 2 )^. 
In this way, equation (7) simplifies to 

U a (n l +an 3 )+n 2 |f = ^~( n *+ aft 3)+^=° ( 8 ) 

where n is the normal to the curve bounding a cross section 
in a plane normal to the x axis. 

All of the subsequent analysis proceeds from Green's 
theorem which relates a volume integral over a region V to a 
surface integral over the surface S enclosing 1 . Green's 
theorem can be expressed in many ways; here it is found con- 
venient to use the forms associated with the linear differential 
equation obtained by equating the left-hand member of 
equation (1) to zero. This results in two different forms of 
Green’s theorem, one for M tt ^ 1, and the other for M m ^ 1 
and prompts the introduction of the following abbreviations 

k=MJ (9) 

^ oo 

If the undisturbed flow at infinity is subsonic (i. e., 
M m ^ 1), equation (1) can be rewritten as 

. c) / ip x c 

+ 22 KtPxVxX l 9 


Figure 1. — View of wing and coordinate system. 


( 10 ) 
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and the corresponding expression of Green’s theorem is (see, 
e. g., ref. 21) 

fjjwm-mmv — <»> 

» 2 


where S2 and <// are arbitrary functions and L(il) is defined as 
follows 

L(n)=p 2 Q zx +U n +Q zl ( 12 ) 

and d9. dv is a derivative along the conormal and is defined by 


dfi 

dr 


dq 

dx 


0 2 >h + 


dn 

by 


n *+ dz n 3 


(13) 


where n 2 , n 3 are the direction cosines of the normal to the 
surface drawn into the region V. 

If the undisturbed flow at infinity is supersonic (i. e., 
M a ^ 1), equation (1) can be written as 


~ fi 2( Pxi + <Pu V + <Pz 2 = k(p x <p xx = k ~ ( 1 4) 

and the corresponding expression of Green’s theorem is 

JJ(*g-»g)<s o« 

r 2 

where the following definitions 


and 


apply. 


L(D)^-/? 2 tt rz +^+a 22 


512 

~bv 


d!2. 2 , 5S2 

d.r 0 ~" l + c)y' ,2+ c)z n * 


(16) 


(17) 


DERIVATION OF INTEGRAL EQUATIONS FOR 
TRANSONIC FLOW 

In this section, integral equations corresponding to the 
transonic differential equation are derived for subsonic and 
supersonic free-stream conditions. One of the principal con- 
tributions here evolves from the attention given to the shock 
waves, or discontinuity surfaces, appearing in the flow fields. 
It will appear (see eqs. (23) and (30)) that the perturbation 
velocity potential can he expressed, for M m less than or 
greater than 1, as the sum of integrals that show no explicit 
contribution of the shocks. ( ’loser analysis of these integrals 
reveals, however, that discontinuities in velocity can appear 
and that they automatically satisfy the shock-polar relations 
(see eqs. (34) and (37)). This section is prefatory to the 
formulation of the transonic integral equations for the par- 
ticular cases of a slender body of revolution and a thin wing. 

INTEGRAL EQUATION, A/co g 1 

The function \p in Green’s theorem, equation (11), is now 
identified with the fundamental solution 1/a of the differen- 
tial equation L(\//) = 0 and the function y is replaced by 12, 
the perturbation velocity potential of the flow field under 


consideration. From equations (10) and (11), the following 
relations hold 



2 


-JJT 


L{<p)dV 




(18) 


where 

a-[(x--x i y+P 2 (y--y l ) 2 +p 2 (z--z l ) 2 ]M (19) 


In these equations the running coordinates in the integra- 
tions are y Y , and z x and <p is to be calculated at a point P 
with coordinates x, y, and z. 

Equation (18) is now applied to the infinite region V 
surrounding the given object to be studied. Some care 
must be exercised, however, in fixing the enclosing surface 
2 since Green’s theorem requires that singularities and 
regions of discontinuity must be excluded from the domain 
of integration. It is to be noted, first, that a vanishes at 
x=x u y=y i, and z=Zy and the effect of the resultant sin- 
gularity can be determined only after the field point is 
enclosed by a neighboring surface and the region V taken 
external to this surface. Second, since shock waves are 
to be expected within the flow field and discontinuities in 
the perturbation velocity components occur across these 
waves, the boundary of Fmust also be drawn so as to exclude 
such discontinuity surfaces. 

In figure 2, a schematic indication of the body and the 
region of integration is shown. The complete three-dimen- 
sional extent of the body has not been pictured; it suffices, 
however, to state that the surface 2 (shown dashed) is 
composed of a sphere of large radius which forms the external 
boundary of V } a sphere of infinitesimal radius surrounding 
the field point P, and a final surface enveloping the object, 
its wake, and its shock waves. 
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If equation (18) is applied to this region and the a priori 
assumption is made that the perturbation field attenuates 
sufficiently fast with distance to negate the contribution 
of the surface integral over the large sphere in the limit as 
the radius goes to infinity, the following expression results: 


?(**,*)=— jj (' G 

0 ~\~ IV 

> -sSSQ r-4 ; ir, T dv <*> 


d 


In tills equation, the integration region over the surface of 
the object and its wake is denoted by 0-\-W . The deriva- 
tives in the surface integrals are, in all cases, along lines 
directed away from the integration surface and into t lie 
three-dimensional domain V since, as follows from equation 
(13), the direction numbers have the same sign as the direc- 
tion numbers of the true normal and the Mach number 
effect is limited to a foreshortening of the longitudinal 
dimension. On the shock surface X the conormals are 
directly opposed on the upstream and downstream faces 
X a and X ft . On the body itself, the conormal derivative can 
be simplified in the manner used in developing the boundary 
conditions of equations (8); that is, from the restrictions 
imposed on the gradients along the body surface, the direc- 
tion of the conormal becomes effectively that of the normal 
n lying in the crossplanes x Y = const. Thus 


d d d __ d 

dv n% &y' n *'bz dn 


on the surface of the body and wake. 

If the triple (spatial) integral of equation (20) is integrated 
by parts r-wise, the resultant form is 



v 


Equation (21) is of particular interest because the integrals 
over the shock surfaces may be shown to vanish. In order 
to prove this, one notes first that the two integrals extend 
over the same geometric surface but that the integrands are 
evaluated, respectively, on the upstream and downstream 
faces and, by definition, the directions of the conormals are 
opposed. When the integrands are combined, the total 
integrand can be expressed as the difference of the two 
terms, one of which contains the factor (^>)x 6 — (<p)\ a and the 
other contains the factor 


^x 1 ^ 2 ^1+^J/ 1 ^2 + ^ 1 ^3“2 < t >X * n * S ) X 

(^P Xl P 2 ni^(PyJl2-\-(PzJlz 2 < f >x \ U l^) x 

The first of these factors vanishes by virtue of the fact 
noted previously that the perturbation potential is contin- 
uous across a shock surface. The second factor can be 
rewritten in the form 

/ 3 2 (ttx 4 -ttxJ(»,)x 6 + K-«\)(n 2 )x 6 + («\— WxJ(« 3 )x 6 — 

X K 2_w \. 2 )( n ‘K 

Also, the change in the velocity vector occurring at the 
shock surface must be in a direction normal to the shock. 
This implies the relations 

(ni) V (n 2 )v (n 3 \= (ux-ux J : (flx -i\) : 

Thus, the second factor to be evaluated becomes 

0 2 («x t -WxJ 2 + K-0X„) 2 + (Ux b 2 -U K 2 ) (ux-u K ) 

[(«**-«0 2 + K~ w O*+ ( w \~ w k)T 2 

The numerator of this fraction, however, corresponds to the 
shock-polar conditions of equation (3) and the expression 
vanishes. 

It finally remains to remark that in the surface integral 
over the body itself, the term ni<p x 2 resulting from the 
integration by parts is of higher order than the normal 
derivative of the perturbation potential. The term can 
therefore be neglected and equation (21) becomes 


<p(x,y,z)=- 



0 + w 



Equation (22) provides another integral expression for the 
perturbation velocity potential, for 1, in transonic 

flow theory. The first integral on the right is algebraically 
equivalent to the expression for <p(r, y, z) in linearized theory 
and the spatial integral is a contribution brought about by 
the nonlinear term of the basic differential equation. It is 
of interest to remark that a derivation ignoring the existence 
of the shock waves can also lead to the same form of the 
equation. In this respect the relation is not unlike cases 
arising in linearized supersonic theory where it becomes 
necessary to study the contribution provided by the fore- 
most shock wave induced by the body. For the majority 
of cases of practical interest, it can be shown that compen- 
sating terms arise and that the discontinuity surfaces are 
taken care of by a formal development that ignores the 
existence of these surfaces (see, e. g., refs. 22 and 23). 
It is not possible, however, to ignore so completely the 
existence of the discontinuity surface in transonic flow and, 
as will be seen in the later developments, equation (20) is, 
for certain purposes, preferable to equation (22). 


428347—57 2 
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Equations (20) and (22) will now be written, for purposes 
of reference, in the following final forms 


V 

(23a) 

£JJJ 

V 

(23 b) 

-•><*,!/, 2>-j; 


im:b iv 

V 

(23c) 


(23d) 


In each of the above equations, <p L has the analytic repre- 
sentation it has in linear theory. The first two relations are 
obvious repetitions of equation (22). The two latter rela- 
tions are transcriptions of equation (20) where the notation 
A(d<p/dv) = (d<p/dv)\ a -\-(d(p/dv)\ b has been introduced, the 
continuity of at the shock surface has been used, and 
where in equation (23d) the variable 


u(x,y,z)=u L (x,y,z)+^r 2 

P " 


_k 

47 r 


JM&> 


(24) 


A detailed account of the application of the two-dimensional 
form of this equation to the calculation of airfoil pressure 
distributions has been given by Spreiter and Alksne (ref. 15). 


INTEGRAL EQUATION, Mco^l 

Use is now made of Green's theorem as expressed in equa- 
tion (15) and 12 is set equal to the perturbation velocity 
potential (p(x,y,z). The direct analogue of t lie derivation in 
the previous section would require that \p be replaced by 
[(x-Ti) 2 — 0\y— ?/i) 2 — P 2 {z— ztf\ A but this leads to the im- 
mediate introduction of a finite-part technique in the inte- 
gration. For the initial stages of the analysis, \p will be 
identified with the fundamental solution c3 of L(\//) = 0 used by 
Vol terra (ref. 24). From equations (14) and (15), the follow- 
ing relations hold : 



where 


JJj SLMdV 

< 25 > 


co = cosh 1 


x — 

&{(y—yiY+{z—z l )'Y 2 


" sinh_1 y.) 2 +(i— 

(g-*i) ln l*-*il + { ( j-*,) 2 +/3 2 [(?/--2/.) 2 +(z-zi) 2 1 } 
|3—*i| /3[(y — yx) 2 +(z— z Y ) 2 ]' A 

is employed to express the integral equation in a form that 
will be of value in establishing a reduction to the case of 
sonic flow. 

The longitudinal perturbation velocity is given by the 
z-wise derivative of any of equations (23). Consider, for 
example, equation (23a). After first isolating the singularity 
at the field point by introducing the limits x { =x±e, one has 


u(x,y,z ) 

c) , 

= dx‘ pL(X 


*’ s) +'!2 l £ If**' [J l ;)*■ 
I + i (Z , ;) *>]=*<w>+ 

.. k C C 1 2 ( v tdyidzi 

2 it JJ 2 ' Pz ' ( * ;2/l,2l) [e 2 +PKy—yiY+P 2 (z—ziy] 3/i ~ 


k 

47 r 


In the limit as e— >0 the influence function in the integrand of 
the double integral is effectively a two-dimensional pulse 
function at the point yi=y, z x =z and of strength 2t/P 2 . 
r rhe expression for u then becomes 


The successful application of equation (25) to the transonic 
problem_hinges on the proper choice of the three-dimensional 
region V and its enclosing surface 2. Discontinuities in the 
velocity components are again to be taken into consideration 
at the shock waves. Furthermore, the fundamental solu- 
tion of Vol terra becomes infinite at yi=y, z Y = z, that is, 
everywhere along the line passing through the field point P 
and parallel to the x axis. Figure 3 indicates the disturbance 
field of the object as well as V and 2 (shown with dashed 
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lines). The bow shock fixes the foremost extent to the dis- 
turbance field and 2 lies adjacent to it and other possible 
shock surfaces as well as the surface of the _object and its 
wake. The downstream limits of the region V are fixed by 
the forecone with vertex at P and determined explicitly by 
the relation 


<p(x 


a r JJ (“I’.-'v)''- 


+ 


r 


* d 

2tt d.r. 


JB--(S) 


dV 


V 


(29) 


(x-z l )=p[(y-y l y+(z-z l ) 2 P (26) 

The inner boundary of \ ' is the cylindrical surface of infinites- 
imal radius given by the relation 

(y~y i) 2 +(z— £ i ) 2 =€ 2 


Equations (27) and (29) may now be written in the various 
forms 

<p(x,y,z) =<PL(x,y,z)+^ ~ jjj 9 ^ 1 2 ( 30a ) 

V 


The conormal derivative is defined by equation (17); on the 
infinitesimal cylinder its direction is parallel to that of the 
normal to the surface, and on the forecone from P the 
conormal is directed along the surface itself. Formal 
analysis yields the expression 



V 


where integrals over the surface of the body and wake are 
denoted by r, over the two sides of t lie sjiock surfaces by \ a and 
A & , and over the enclosed volume by V. In each case, only 
that portion of the surface or volume lying within the forecone 
of P is included in the integrals. The surface integrals over 
the forecone itself vanish because aJ and dco/57 are zero. It 
should be noted that the forecone is that of linearized theory 
and has no rehi ionship to the region of dependence in the 
actual flow field. 

Integration by parts, in the last integral, leads to the 
relation 


*<*,!/,*) 2„S/J ["(™ _ 2»'>) _ »li"} i - _ 

is (If + J] 


+ 


(28) 


V 


Equation (28) is the form, for ^ 1, analogous to equation 
(21), for ^ 1, and on the body and wake surfaces involves 
the approximation 


d d , d d 
d7~^ 2 c)?/ dn 


It is not difficult to show, from the shock-wave relation of 
equation (3), that the combined integrals over the surfaces 
A a and A & vanish. The perturbation velocity potential can, 
therefore, be given alternatively as 


V 

-«(*«, Z)-^ 4 (|f)‘ as - 

V 

-«(W)-^//a(£?) = <S- 

X 

MTiY d 
m 2rJJJ \dXi 2 )a 


(30b) 


(30c) 


dV (30d) 


V 


where <p L has the same analytic form as in linearized super- 
sonic theory and use has been made of the relations 


o = cosh 1 


X — Xi 

0[(^-2/i) 2 +(2-2l) 2 P 

—1 


dco _ 

d*i “ [(x—x l ) 2 —P\y—yx) 2 —P\z — z x ) 2 } * “ 


(31) 


Comparison of equations (23) and (30) shows once more 
a difficulty that appears in linearized analysis of subsonic 
and supersonic flow, namely, that complete parallelism be- 
tween the formulas is not achieved directly. This parallelism 
can only be established after interchanging the order of 
integration and differentiation and, because of the singulari- 
ties involved, it becomes necessary to introduce the concept 
of finite-part integration. Furthermore, it is well known 
that the resulting multiple integrals can no longer be written 
in a unique form (see, e. g., ref. 25) but must be expressed 
differently, depending on t lie order in which the integrations 
are to be performed. No attempt will be made to develop 
these ideas further at the present time. 


CALCULATION OF CONDITIONS ON SHOCK SURFACES 

It is of some interest to study equations (23a), (23b), (30a), 
and (30b) as the field point approaches a discontinuity sur- 
face and to discover the mechanism by means of which these 
basic equations furnish the velocity jumps associated with 
the shock waves in the field. To this end, consider first the 
case M m ^ 1. Figure 4 shows the geometry of the problem. 
The bow wave induced by an aerodynamic shape is indicated 
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and the point P, at which conditions are to be calculated, is 
chosen arbitrarily close to the rear surface of the wave. The 
surface of the wave can be replaced locally by a planar (‘le- 
nient and a new coordinate system £, 77 , f introduced with 
the origin fixed at the intersection of the line 1/1 = 1 /, Zi = z, 
and the bow wave. Point P then lias t lie coordinates £, 0, 
0 and the planar surface is given by the linear relation 


+ + — 0 (32) 

Since the bow wave is situated upstream of the linearized 
disturbance field, u L is zero and, from equation (30a), the 
perturbation velocity is 

*<*■ “• °>-S § JIf T (M) * <*» 

where aj = cosh” 1 [(^—^i)/jS(7/ 1 2 +i'i 2 )^]. By virtue of the field 
point’s nearness to the bow wave, the term u 2 /2 in the inte- 
grand is assumed a constant and one then gets 


where 


n n\ k u * d 2 r Z2 , . r Y * j r X2 da 

0. ' »)-2i T 5? J * l,s ‘ J „ d " J .V, i W, di ' 


— — (&»7l+cfi)/<l, X 2 — £ — + 

xr v _&(a£+cfi)±a/3[a 2 £ 2 +2ac£fi— (a 2 /3 2 — 6 2 — c 2 )£i 2 p 

~ ‘ „2o2_7.2 


z„z 2 = 


a 2 |8 2 — 6 2 

_fl|[c ±(ffl 2 |8 2 — 6 2 )^] 
a 2 /3 2 — 6 2 — c 2 


Integration and differentiation yields 

/J- A A\ O' 

( *’°' 0) ~ 2 — a 2 0 2 +6 2 +c 2 


( 34 ) 


It remains to show that this result, derived from the 
integral expression for the perturbation velocity potential, 
is consistent with the result one would get from the shock- 
polar relation of equation (3). At the downstream face of 
the bow wave, equation (3) becomes 


— P 2 u b 2 + v b 2 + W = ku b 3 /2 


The incremental velocity vector occurring at the shock sur- 
face is, however, normal to the shock surface and this yields 
the relations u & :r 6 :ic & =a:6:c. Substitution into the shock 
polar relation gives 


ku 1 \_k /J o? 

2 \P 2 Ub 2 — r,r~ w, 2 ) 2 —a 2 (3 2 ■ b 2 c 2 


(35) 


in agreement with equation (34). It therefore follows that 
the integral expression, for 1, will adapt itself on the 

shock surface to any bow wave consistent with given body 
geometry. This result can also be extended to include any 
shock wave in the flow field. 

An analogous procedure follows for the case M m ^ 1. Let 
the shock surface in the vicinity of the point P at £,0,0 be 
given by equation (32) and assume that u 2 is composed of a 
continuous part and a discontinuous part that has the con- 
stant value u 2 ahead of and u b 2 behind the shock. Equa- 
tion (23a) then yields 


lim [u{£ — ,0,0) — ^(£+,0,0)] = 
£-0 


lim •— ( u<?—u b 2 

£->() 47T 



ad^dry, 

M+b Vl +ct l ) 2 +a 2 t3W+a 2 p 2 tf] l/2 


(36) 


where the double integration extends over the region of 
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discontinuity. If the differentiation with respect to £ is 
now carried within the integral signs and £ allowed to ap- 
proach zero, the value of tiie integral becomes independent 
of the original limits of integration. In this way one gets 


u a —u b =lim [u(£— ,0,0) —^(£+,0,0)] = 


ka 2 


u 2 —u b 2 




a 2 0 2 +6 2 +c 2 


(37) 


It can be shown, as previously, that equation (37) agrees 
with the result given by equation (3), the shock polar con- 
dition, for ^ 1. 


REDUCTION TO SONIC FLOW THEORY 

In this section, the previously determined equations will 
be studied in the limit as sonic, free-stream speed is reached. 
The integral relations then assume forms that correspond to 
the nonlinear differential equation when 0=0. 


INTEGRATED STRENGTHS OF EXTERNAL SOURCES 

It is proposed here to determine a relation that will prove 
useful in the following section in connection with the reduc- 
tion of the integral equations to the special forms appropriate 
for M m = 1 . This relation will be recognized subsequently 
as connecting the integrated strengths of the exterior correc- 
tive sources in the crossplane x=x 0 = const, and the rate of 
change of body cross-section area. As a means to this end, 
equation (10) is written in the form 

t (38) 

where V 2 is the two-dimensional Laplacian operator in the 
transverse plane. Each term is then integrated over the 
entire x=x 0 plane external to the body. The double integral 
involving V 2 <p can be partially integrated and converted 
into a line integral by application of Green’s theorem for a 
plane 

(©„,*• m 

where the line or curvilinear integral extends around the 
curves C enclosing the region of integration of the double 
integral. For cases in which the plane x=x 0 does not inter- 
sect any shock waves, the region of integration can be taken 
at once as the entire x=x Q plane exterior to the body. If the 
assumption is made, as in linear subsonic theory, that the 
normal gradient of <p attenuates with lateral distance suffi- 
ciently fast to suppress the contribution of the curvilinear 
integral along the outer boundary, the boundary conditions 
of equation (8) permit one to equate the line integral along 
C to —U a> S / (x 0 ) where S'(x 0 ) denotes the longitudinal 
gradient of body cross-section area. The integrated form of 
equation (38) thus becomes 


Jiu tl., dvd!=$ iKsi), 


dydz—U„S'{x 0 ) (40) 



Figure o. — View illustrating intersection of shock wave and 

x=x Q plane. 


ties occur which require that the integration region must be 
divided into two parts, one lying between the body and the 
shock wave and the other extending beyond the shock wave 
to infinity. Application of Green’s theorem to each region 
and addition of the separate contributions results in addi- 
tional line integrals carried around the two sides of the shock 
surface. These two line integrals can be combined into a 
single line integral, in which case equation (40) can be written 
as follows 


k J.(Qi CL.^+l [ A <ML 

F [J ^ dydz—U„S'(x 0 ) (41) 


where A(dv/c)n) = i><pal()n — ()(p b ldn and where in the single 
line integral the integration extends around the curve de- 
scribed by the intersection of the shock wave and the x=x 0 
plane, and the normal n is taken as directed away from the 
body. 

In the subsequent work, attention is to be directed toward 
results at M m — 1. We assume here that in the limit as 0— >0 
the first term in the right-hand member of equation (41) will 
vanish and one then gets 


*JJ(5tL** + J'( a §S). 


dax=-U„S'(x 0 ) (42) 


It will become evident in the discussion contained in the fol- 
lowing section that the left side of equation (42) represents 
the integrated strengths of the exterior corrective sources in 
the cross-plane x=x 0 = const, in the limit as M m approaches 
1. It follows that under conditions corresponding to sonic 
flight speed, the total source strength in any transverse plane 
is zero: The sum of the sources within the body or wing 
(sources appearing in the term <p L ) is of equal magnitude but 
opposite sign to the corrective sources required by the non- 
linear term in the differential equation. 

Equation (42) allows one to make some conclusions about 
the lateral attenuation of the r-wise gradient of u 2 . Consider 
the double integral as written in terms of polar coordinates 



In the more general case, however, in which the plane x=x 0 
intersects a shock wave as illustrated in figure 5, discon tinui- 
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Since the definite integral must converge, it follows that the 
integrand attenuates faster than 1 /r x 2 and if the assumption 
is made that the integrand has a purely algebraic character, 
one concludes that for large i\ 


0 U 2 _J_ 

bx 2 ~~n 2 +* 


(43) 


where N is some positive constant. 

A check on equation (43) is provided by the work of 
Guderley and Yoshihara (ref. 26) on axially symmetric flow 
at sonic speed. In that analysis, for large r lf 


0 //-' 1 
dx 2 7 


and this is in agreement with equation (43) when A r =2/7. 
r Phe same reference also gives 

b<p 1 

d/w, 977 


which serves to substantiate the assumption made earlier 


that the curvilinear integral along the outer boundary can 
be neglected in equation (40). 

INTEGRAL EQUATION FOR SLENDER BODIES. M*>=1 

In the reduction of the integral equations to the case £=0, 
methods analogous to those employed in references 25 and 
27 will be used. Attention will be confined here to field 
points at a finite distance from the body so that in the limit 
as fi approaches zero, the term fir can be assumed to approach 
zero uniformly. As in conventional slender-body theory for 
linearized flow, t lie longitudinal distribution of cross-section 
area S(x) is assumed to possess a continuous x-wise deriva- 
tive. The method of reduction can be exhibited in a suffi- 
ciently general form if a lifting body of revolution is con- 
sidered. When 1, the perturbation potential for a 

body of revolution <p B follows from equations (23) and for 
the purposes at hand the form (23d) is preferable. For 
sufficiently slender and smooth bodies the term p L can be 
expressed in terms of a rectilinear source and doublet dis- 
tribution and the corrective source distribution appearing in 
the triple integral then extends over all space external to 
the x axis. If k(x ) is the x-wise distribution of doublet 
strength, equation (23d) becomes 


M (rr aw U <° r S\x l )dx l sill 6 d_ f°° (z—xMxJdXj 1 d_ C C ( <V\ (x— xQ , 1 s — a?i | + [(a — x x ) 2 +fi 2 p u 2 ] ** ^ 
4 * Jo [(x—Xi) 2 +fi 2 r 2 ] H 4tt r dxj 0 [(x—x x ) 2 +fi 2 r 2 ] y * 47rdxJJ\ dv) B \x— z x \ " fip n 


k d CCC ( 5 <Px *\ (x— xQ 
4 t t d.rJJJ \dx l 2 ) B \x—Zi\ 


) ln \x—z x \ + [(a— x 1 ) 2 +fi 2 p Il 2 }y* 


fip i 


rxdridxiddi (44) 


In equation (44), cylindrical coordinates x, r, 6 are used where 
r 2 =y 2 -\-z 2 > 0=tan _1 ( zjy ) and the notation 

Pn = [r 2 + r{ 2 — 2rr l cos (0 — Q x )]X 


has been introduced. 

Under the imposed conditions, the first two terms in the 
right member of equation (44) reduce, as shown for example 
in reference 27, to approximately 


-S\x) Inr- 


k(x) sin 0 U 


Zrr 


4t 




X—Xi 
X — Xi 


j 2\x—Xi\j 

in - — — dXi 


It remains to attempt a corresponding modification of the 
two remaining integrals. Consider, next, the triple integral. 
Since the integration extends over all points in space, one 
encounters a nonuniformity of convergence in the logarithmic 
influence function when fi becomes vanishingly small and r\ 
becomes infinitely large. From equation (43), however, it is 
known that du 2 /dx attenuates rapidly with increasing r l and 
the resultant error in miscalculating the effect of the loga- 
rithm for large r x is thereby reduced. The triple integral is 
therefore approximated by 


— — f f f u ^ x uruB x T \ x ~ x \ 

2 tt dxJJ J Ld-Pi 2 J \x— x x 


In 


2\x—Xi\ 

fipn 


dx l r\dr ] dB\ 


This term can be rewritten as 



u B \x y r u e x ) 


J In pnridriddi— 


k b C x— x x , 2\x—xi 
47 r dxj |x — X\ | n fi 



/ U' B 2 (xi,ri,6i) 


J n dr v 


dO x 


Consider, finally, the integral over the shock surface. For 
fi near zero, one gets the expression 


-LUK 




i n 2 \ x - x ' \az 

)\X-Xx\ ln Pru dx 


which becomes 


sj> I? 


( x >yi> z i)lnpnd<T \ — 


1 d r x—x x 
4 t dxj \x— x x \ 



(x,y l ,z 1 )dax 


where da\ represents an element of arc on the curve deter- 
mined by the intersection of the shock surface and a plane 
normal to the x axis. The normal n, as in equation (42), is 
directed away from the body. 

The reduced form of equation (44), for the case M m = 1 , 
is given by the sum of the above three expressions. In the 
forms given, an apparent dependence on fi remains in the 
expression 


1 d C x—x x 
47 r dx J \x—x { \ 


, w 2jy^,[^, (a . )+ 



2 


r 1 dr 1 d0 l + 



d<p B 

dn 
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The bracketed term vanishes by virtue of equation (42), 
however, and the dependence on /3 disappears. The sonic 
form of t lie integral equation for the lifting body of revolu- 
tion is, therefore, 


, U^S'ix) 
<p B (x,r,6)= — Inr- 

Z.7T 


k{x) sin 6 


w7T T 


27 r j\ 


A + lnpnfa+ 
X on 




(45) 


Starting under the assumption that and using 

equation (30d), one can derive the same result, the only 
essential difference arising from the fact that in the limiting 
process care must be taken to restrict the disturbance region 
V to that portion of space within the Mach forecone from 
the field point at x,y,z. 

Equation (45) expresses the sonic equation in the form 


<PB(x,r,6)=<p2 B (x;r,e)+^ f. A lnp n d<r x + 


£//(*¥) lnpnridriddi 


(46) 


where V 2 u (x)r,Q) is the harmonic potential for the body of 
revolution in transverse planes. As expressed, this sonic 
form of the integral equation for the perturbation potential 
is identical to the integral equation corresponding to the 
transonic differential equation at p=0] that is, it gives, for 
a flow field in which shock waves may possibly occur, the 
integral equation corresponding to the partial differential 
equation 


i 7+1 

<P Vl)\ i PzZ JJ tpxtyxx 


(47a) 


and admits discontinuity surfaces for which the difference 
relation 

7 — f - 1 iDr — I - (pr 

(.<Pv a — Vvy-\-{<Pz — <%) 2= 77 ^2 - <Px b ) 2 (47b) 

is satisfied. The direct derivation of equation (46) would 
follow from an application of Green's theorem in the trans- 
verse variables to equation (47a) without the introduction 
of the slender-body assumptions. 

The interpretation of equation (42) in terms of net source 
strengths is now apparent. Each of the three terms in the 
right member of equation (46) provides two-dimensional 
sources in each transverse plane ar=const. and the perturba- 
tion flow field is simulated by the combined effect of these 
sources and in lifting cases, a doublet term. (The doublet 
term is of no concern in the present discussion since its net 
source strength is zero.) The first term, contains a 
source on the x axis with strength fixed by the gradient of 
area, S'(x); the second term represents a curvilinear distri- 
bution of sources around the shock wave with strengths 
fixed by A(d<p/dn) ; the last term represents a planar distribu- 
tion of sources with strengths determined by the nonlinear 
term (k/2)(d/dx)<p x 2 . Equation (42) thus states that at 
M a> = l the combined source strength must vanish in each 
transverse plane. 

Equation (46) corresponds, for the body of revolution at 
M m = 1, to equations (23c) and (30d) in that it contains an 
explicit contribution from the shock wave and is expressed 
in terms of the basic singularities of the differential equation. 
A form analogous to equations (23b) and (30b) can also be 
derived as follows. Let r=R(x,0) be the equation of the 
shock surface. The relation 


A f 

2wJ 0 


~ U "- Inpnr^dffi- 


R(x,di) 






dd i 


then holds. Substitution into equation (46) yields 


<p B (x,r,d)=<p 2 B (x;r,d)-^ - 


dlik 2 
u B 2 


0<Pb_ Ox 2 

0/z 




/«pii<7<7 X + 9 ^JJ' lfLpii/‘\dr \d6\ 


and through use of equation (47b) one gets 

Pb(x,v fi)=(p2 B {x]r , 0 )+ 2 ^ dxJ'J 2 lnpnridriddi 


(48) 


The above results have been worked out in some detail 
for the body of revolution. The sonic equations for other 
shapes follow similarly through a reduction of the general 
transonic integral equations or can be expressed directly 
through consideration of the sonic differential equation. 
The final equations for both cases appear as follows: 


k d C C u 2 

<p(x;y,z)=<p 2 (x;y,z)+ 2 ~^ I ~~ lnp u rxdr x d& x (49a) 
=<pi(x; y,z )+ ^ A lnp r ,d<Tx+ 

2tt fji dj ~2~ I (49b) 
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where <p 2 (x;y,z ) is a two-dimensional harmonic function 
which, for the body of revolution and for the planar wing is, 
respectively, 


U m S'(x) j k(x) sin 0 

<P 2 b = — ^ Inr — ^ 

B lir 2 wr 

I rs 2 w 

Aiv w (x;y l )ln\(y~y l ) 2 +Z 2 } y2 dy l + 

Z7r J -S l (x) 

1 P 2 '*’ A / \ zdyi 

^J^ lM Avw{x ’ yi) [(y-yir+z 2 ] 


(50a) 


(50b) 


In tlie latter expression, 


Aw(x; yO=w(x; yi) 2 ^—w(x; yi) 2=C - 

A<p(x; y x )= <p(x; yi) 2=0 +— <p(x; yO.-o- 

and the lateral boundaries of the wing plan form are fixed by 
— Si(x) and s 2 (x). 


SLENDER-WING THEORY IN LINEARIZED FLOW 

In the preceding development, t he integral realation for 
the perturbation potential in sonic flow has been expressed 
in a form that follows from an application to equation (47a) 
of Green’s theorem in transverse planes. The determination 
of a solution thus depends to a large extent on the evaluation 
of the effect of the two-dimensional singularities that are 
placed throughout the exterior portion of the flow field. 
Examples of a direct attack on a similar problem are to be 
found in the calculations of two-dimensional transonic flows 
by Oswatitsch, Gullstrand, and Spreiter and Alksne, refer- 
ences 9 through 15. In the present report, an indirect 
attack is to be made, following the ideas of Whitcomb and 
Oswatitsch (refs. 6 and 7) by relating the solution for a 
slender wing to that of a body of revolution. The analysis 
will show that once one establishes the details of the crossflow 
potential fields associated with a wing and its related body 
of revolution, the residual disturbance fields near the two 
bodies are the same to a certain order of accuracy in terms 
of the slenderness ratio. The mechanics of such an approach 
can, in fact, be observed in linearized wing theory and such a 
development will be given in this section as a prelude to the 
subsequent sonic theory. Attention will be limited to the 
subsonic case and, as an added simplification in the analysis, 
only wings possessing lateral symmetry will be considered 
although such a restriction is not essential. 


ANALYSIS 

The linearized equation for subsonic potential flow is 


^P UU ^P ZZ (51) 

and if Green’s theorem is applied, formal manipulation leads 
to the following integral relation for the potential <p w of a 
planar wing 

<Pw(x,y,z)=<p iw (x;y,z)—~ ij j ^^(x.r^OM^+r, 2 — 

2rr x cos (0— 0 1 )PVitfr 1 (/0 1 (52) 


No integrals along possible discontinuity surfaces are neces- 
sary since shock waves do not appear in linearized subsonic 
flow theory. Equation (52) is linear in y and can be sep- 
arated into additive expressions contributing to <pw % %, the 
potential associated with the thickness distribution, and 
<Pw,a, the potential associated with camber and angle of 
attack. In this way one gets for the perturbation velocity 
components 


u w ,,=u hl . <Pxx w Jn[r 2 +r i 2 —2rr l cos (0— 0,)PV,oM0 1 


(53a) 


Ww,t=W-i w — I* <Pxx w Jn[r s 4-ri 2 — 2rr x cos (0— 0,)pr 1 c?r,rf0, 


and 
u 


(53b) 


8‘ 2 5 rr 

w,a=u 2w — J <Pxx w Jn[r 2 +r i 2 —2rr l cos (0— 0i)]*V l c6yJ&i 

(54a) 

• §§Vxz w ln[r 2 +ri 2 --2rr l cos (0— 0 1 )]^ 1 cM0 1 


(54b) 


The corrective integrals in equations (53) and (54) ob- 
viously do not modify the area distribution in the thickness 
case by virtue of the vertical symmetry in the flow field nor 
the load distribution in the camber case by virtue of the 
vertical asymmetry in the field. It follows that 


AW\y i Aw 2}v t , A r u^ ta — Au 2w a (o5) 

where the delta notation denotes the increment in the func- 
tion in passing through the plane z=0, that is, the difference 
between the values on the upper and lower surfaces of the 
wing. As a consequence, <p 2w can be expressed in the form 


<P2 w (x'>y> z ) 

u a 


-if. 


s(z) AWi 


w, t 


S{x) 


ln(r 2 +yi 2 — 2ry x cos 6)^dyi+ 


J_ f s(z) A<p Wta r sin 0dy x 
2? rj— , (i) U m r 2 +y\— 2ry x cos 0 

If exact conditions on the wing surface are to be sought in 
linearized theory, equations (53a) must also be satisfied. 
For example, in t lie direct case of given thickness, equation 
(56) predicts u 2w t and equation (53b) is then used to deter- 
mine the exact streamwise velocity component as affected 
by the external-source integral. In the direct case of given 
loading, equation (56) predicts w 2wa and equation (54b) is 
then used to calculate the true wing camber, modification of 
w 2wa being produced by the integral term. The difficulties 
of such calculations are so disproportionate to those of solv- 
ing the linearized equation by standard methods that they 
appear to add needless complications to a relatively simple 
problem. In transonic theory, however, the right-hand 
member of equation (51) is replaced by a nonlinear term and, 
in the absence of more obvious methods of attack, the diffi- 
culties involved in such an approach become less of a deter- 
rent; in particular, the integral forms of the corrective terms 
are of added interest since they are suited to approximations. 
The details of such an approach will not be considered further 
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at the present time since for slender wings the use of a re- 
lated body of revolution yields information of a sufficient 
order of exactness in both linear and nonlinear theory. 

In the following work, a complete knowledge of the solution 
for a body of revolution will be assumed known. Thus, if 
S(x) is the cross-sectional area of the body and the distribution 
of lift is fixed by k(x), doublet strength per unit of length, 
the linearized solution for the body of revolution is 

<p B (x,r y $)_ 1 r S\x x )dx i 

U m 4tt Jo 

P 2 r sin 9 K(xi)dx x , ~ 

47T Jo 

and the perturbation potential in the transverse plane is 


r> 9) S\x)ln r k(x) sin 0 

U a 27 T -7T/- 


(58) 


The integral relations for the wing and the body combine 
to give 


<fw 


-<Pb=<p* w —V 2 b ~- 2 t jj (<Pw—<pb)zx ln[r 2 -\-r 2 — 

2rr x cos (6— 9ii\ X/ ^Vidr x d9i (59) 


Equation (59) is exact, subject only to the restrictions of 
first-order perturbation theory, but with added restrictions 
on the geometry and loading it is possible to show that the 
magnitude of the final integral is negligible to a certain order 
of accuracy. We now assume the wing is slender, that is, 
s(z) is small in comparison with over-all wing length. Let, 
furthermore, the wing and body be of equal length and have 
identical longitudinal distributions of cross-sectional area. 
This implies 



Aw w (x,y x )dy x 

U m 


(60) 


and establishes the condition that the distribution of two- 
dimensional source strength in equation (58) is equal to the 
strengths integrated in the transverse plane of the sources 
appearing in equation (56). It will also be convenient to 
equate in the same manner the doublet strengths in those 
two equations and one is led to the relation 


k(x)= 


J s(z) 
-Six) 


A<p 2w (x,y x )dy x 

U m 


(61) 


The first objective will be to show that for field points in 
the vicinity of the slender wing the first two terms in the 
right member of equation (59) differ from the left member 
by an amount that is of higher order in s/l. The evaluation 
of the error term can be performed by an iterative process 
in which the first step starts with the approximation 


<pw~<p2 w <P2 b ~\~<Pb (62) 

Before integration, the integrand in equation (59) will be 
written as the product of Fourier expansions of its terms. 
For the logarithmic term, one has 


ln[r 2 +r x 2 - 
2rr x cos (0— 0j)]^= -< 


and, similarly, 


Inr i 


Inr 


__ ^ / r\ m cos 
£?i \rj 

771= 1 \ 1 / 


os m{6—9i) 


m 


r^r x 


jos m(6—d 1) 
m 


(63) 


-J T\ ^ 1 


r sin 6 


r 2 Ay x 2 —2ry x cos 6 


|£(*)" sin m«, |»|£. 
y 1 77i=i \ r / 

1 00 / r \ m 

-ZlL-j sin md, r^\ Vl 

y i w=i \y i/ 


(64) 


If equations (63) and (64) are used, together with equation 

(56), 

and conditions of bilateral symmetry are imposed, one 

gets 

<P 2 w(x;y,z ) S\x) ^ a 2m (x) cos 2 md 

U a l 2 rl 1 (r/sr h 



b- 2 m- i(x) sill (2m — 1 )0 ^ 
(r/s) im ~ l 

(65) 

where 



. — s r w Aw 2 {x,yi) /yiY m dyi 

2m() 4 rmlJ. Hx) U w \sj s 

(66a) 

and 

b (x) - i r 

6 *"-*(*HmJ.. U) u m WJ s 

(66b) 


Equation (65) holds true beyond the circle of radius s(x) en- 
closing the transverse section of the wing; within this circle, 
the expression is 


f 2 ir(-e;y,g)_ S'(x) 


UJ 


2*1 + ^i(2n+iy 


ter-*] 


+ 


^ cos2w0^ A / \ 

777=1 Mil 77 =0 


1 


2m—2n—l 


(r \ 2n+1 i 

< 1 \ 

W ' 

\2m-\-2n+\J 


+ 


^\2m-\-2n— l"^"2m— 2n— l)J’ /=,s 

where the coefficients A 2n and B 2n are related to the boundary 
conditions through the expansions 


Aw 2 (x,yi)_ 

TJ 

X lA in {x) /; 

|)“ 

(68a) 

u a> 

77=0 S \ 



^i(x,y i) 

UJ 

(*)(« 
77=0 \C 

:T 

(68b) 


Once the coefficients a 2m , 6 2w _ 1 , A 2n , and B 2n are related in 
magnitude to the geometry of the wing, the size of the integral 
term in equation (59) can be estimated. Since Aw(x,y)/U m 
is proportional to t(x)/l , the wing’s thickness ratio, it follows 
from equations (66) and (68) that the following order esti- 
mates hold 

a 2m {x ) = 0 (ts/l 2 ), A 2n (x ) = 0 (ts/l 2 ) (69) 
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Let, furthermore, A (p(;x,y)IU m be assumed proportional to 
a ( 2)9 Or) where a(x) is a measure of local angle of attack or 
camber; equations (66) and (68) then yield for the remaining 
coefficients 


&2m-i(s)=0(as/Z), B 2n (x)=0(pcsll) (70) 

Equation (67) can also be written in the simplified form 


< tew(x\y,z) S' (x) 

UJ % rl 

S sin (2rn—\)QH m (x, 

■m = 1 \ *9/ 


lns+F 


(*'0 


’ cos 2 md n 

+ -?,“2 ^r G " 



with the order estimates 


F (x, 0 =o (ts/n, G m (x, 0 =o (ta/P), H m (x, 0 = 


(72) 


Evaluation of order of error in equation (62). — Equation 
(59) is now written in the form 

<p w (x,y,z)=<p 2H Xr;y,z)—<p 2H (x;r,6)+<p n (x,r,6) + I(x,r,0) (73) 
where 

£2 /*co /‘2 tt ^2 

I(x,r,6)= — ^— 2 ( ( Pw—<PH)l n [f‘ 2 -\- 

■ZTTj 0 Jo OX 

r*—2rr\ cos (0— 0i)fi*/*ie//v/0i (74) 

and an estimation of the order of magnitude of l(x,r,0) is to 
be made for field points in the vicinity of the wing. The 
approximation of equation (62) is to be used together with the 
given expansions of the two-dimensional perturbation poten- 
tials. It suffices to simplify the analysis and estimate the 
order of the error at /* equal to . 9 . One then gets 


I(x,s,0) . I 2 (x,s,d) 

UJ ~ UJ ^ UJ 


£ JX P'£r («■ ?)+, 5 t « (*.?>.?, *• < 2 "- '>»'»- 0 ' 5 )] [*— 


r\dr\dO\ — 


(75) 


After integration with respect to 0 } one has 

fs (-Y"cos 2me+^±r (-) 

UJ 2 J s ,„=] L 2 m \rj 2m + \\rJ 


4w+2 


f‘ r <?■ r/'A” cc,s 2me r" ( . 

•rj,, SLw <,»* M x --)+w 


sin (2m + 

sin (2m — 1)0 
(2m — l) 2 


1)»] [Qy «.i+r (.,?)] iiw,*, 

m ( s, (7 )1 r,dr, 


From equations (69), (70), and (72), the order of magnitude 
of I(x,s,6)/UJ is given by 


i(x,s,e) 

UJ ' 


: /3 2 0 ( - 


(ji Ins) 


(76) 


The first step in the iteration fixes the maximum value of 
the error incurred in neglecting the integral term of equation 
(74). In the vicinity of the wing, therefore, the perturbation 
potential can be expressed as in equation (62) with an error 
of the order given in equation (76). 

Near the body, a further reduction of the difference <p B —<P2 B 
is possible since the explicit equations for the body of revolu- 
tion are available. This yields the usual form of the slender- 
wing solution. Thus, from equation (62), 


<Pw~(P2 W -hf(x) 


(77) 


where 


Q2 U2tt r 00 

f(x)=lim(<p B —<P 2 ») = “s" <Pxx B (x, rjjnriridr 

r->0 -ttJ o Jo 


idOi 


(78) 


From equation (57), f(x) becomes, for the subsonic case, 


/(*)= 


P 2 u m e> 

Air dx 




Inriridri 


8 t r dxjo \X— Xj\ [_ 

= _(£A (* , (* = 

4t dxj 0 \x— 


xtf+pW] 3 ' 2 

} 


In -\-lnA I dxi 


u. a 

Zl\ P 


(79a) 


It is not difficult to show by a similar analysis taking into 
account the possibility of discontinuities in the flow that 
equation (77) holds also for M m y>l and that f(x) then has 
the form 


f(x) 


2tt dxj o 


(xjn 


2(x — Xi) 

P 


dx i 


(79b) 


SLENDER-WING THEORY IN SONIC FLOW 

The extension of the foregoing derivation to the nonlinear 
case will be given in this section. The iteration procedure 
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designed to discuss the linear problem is effectively the de- 
velopment of an expansion in terms of the slenderness param- 
eter of the wing and with appropriate restrictions on wing 
angle of attack can be applied with little modification. 

ANALYSIS 

Equation (49b), the relation fundamental to the following 
discussion, provides that the sonic expression for <p valid in 
the x— x 0 plane can be considered to be composed of three 
terms; <p 2 , a line integral around the possible intersections of 
the x=x 0 plane and a shock wave, and a surface integral 
over the entire portion of the x=x 0 plane exterior to the body. 
In man}" important cases, simplification occurs because the 
line integral introduces no contribution to the values for p 
in the vicinity of the body. Perhaps the simplest case in 
which this situation develops is that encountered very 
frequently at M m = 1 in which the shock waves are situated 
entirely downstream of the most rearward point. A second 
case in which the line integral introduces no contribution 
occurs when the discontinuity surface is situated in an 
x=Xq plane and is, therefore, essentially a normal shock 
wave. The discontinuities associated with the normal 
shock wave are contained in the contribution of the double 
integral. Since most sonic flows about smooth wings or 
bodies probably fall into one of these two cases, attention 
will be confined in the following discussion to those cases in 
which no contribution results from the line integral. Thus, 
if equation (49b) is written first for a wing and then again for 
a body of revolution, and the latter is subtracted from the 
former, the following relation is obtained: 

<^pp-(iC,I/,2) = <P2 W ~f" *Pb (*T p yO) -f - J(X,T )@) (80) 

where 


J ( x >r,0)= YT j o J o 2^ ( ^ 2_Ui?2) 

In [r 2 -r\ 2 -2rn cos {e-ei)]*r u dr u dBi (81) 


The quadratic nature of the integrand in equation (81), 
together with the additive dependence on thickness and 
camber in the transverse-plane potentials, prompts one to 
simplify the analysis to cases involving a thin wing of given 
thickness but limited to an angle of attack or camber- 
length ratio a. that is small in comparison with the thickness- 
length ratio t/l. In this way, sufficient information is re- 
tained to establish the relationship between the wing and 


body flow fields for the thickness case and at the same time 
determine a linear dependence on angle of attack of the wing 
loading in the vicinity of qj= 0. Under these conditions it 
will be possible to relate the wing flow field to that of a. body 
having the same area distribution but not inclined to the 
free stream. Consistent with these conditions, the following 
equations apply: 

<p2 |r =(fV ( +¥V„, ‘Pl,=‘Pl nt l82) 

and the perturbation potentials for the wing and body can 
be expressed as 

<PW=<PW, aj <PB = <pB,t (° 3 ) 

where the subscripts t and a identify the contributions at- 
tributable to thickness and camber. The term Uw 2 u b~ in 
the integrand of equation (81) can now be approximated by 

Uw 2 -u B 2 ~(u 2w -u 2B y+2u B .,(u 2w -u 2Bt ) + 

U * B J U *W.a + 2 U B- <84 ) 

where the initial assumption 

<PW ~<P 2 W ~ <P' 2 b * + <£/*, t (35) 

has been made and higher order dependence on a has been 
deleted. It remains to show, through the evaluation of 
J(x,r,e) y that the assumption made in equation (85) holds. 
It should be noted that, to the order of exactness of this equa- 
tion, normal shocks on the wing and body are situated at the 
same longitudinal station. 

The first two terms in the right member of equation (84) 
depend solely on the thickness distribution of the wing and 
body, and the two remaining terms contain the effects attrib- 
utable to the lift and thickness combined. Substituting 
from equation (84) into equation (81), we see that the first 
two terms contribute to J(x,rfl) a function that is symmetric 
about 3=0 and the two latter terms contribute an asym- 
metric quantity. From equation (80) one then gets 

t = ^^' 2 w t i Allw, a = ^'U' 2 w a 

These relations are identical to those given in equation (55) 
and, as a consequence, equation (56) necessarily remains 
valid along with the expansions given in equations (65) and 
(67). 

We now approximate J(x,r,8) at r=s: 


J(x,S,d)^ J\(x,8,0) 
UJ ~ UJ 


JziXtSid) k d 1*°° f 27r cos 2md\~Y A 0 Tn cos 2m0i . d2m( x ) cos_2? 

-tcr = c5ij, .), \LS — + U=' I w #- 


sin (2v— 1)0, 

(r./s) 21 '" 1 


2 u B , «1S 

s «($)]+»*& l < + r (*- ?) + s W ?)H^+ 


dt'idOi (86) 
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The orthogonality of the trigonometric terms, together 
with the fact that u B<t has no dependence on 0, permits one to 


achieve considerable simplification after the 0i integration. 
The resulting expression is 


J(x,s,0) d 

U / 4 7 r d.r 


( X m(%)~ / v w o ,, r 2in0 (i 2 

J \ '£ (F, £ L nsr- ft 




+ 


sin (2m — 1)0 b' Jm -i(z) 


l! + ' ' }'• 


dr 1 + 


k d 
4 ir d.r 


j s ym 


+ 


(2m — 1) (ri/s) 4 "' - ^ 

S '"2;"'*' ( *■ 7 )J '« 8 r ' M < +L 5 *”«•{* s+ r ('- 7 )] ln 8 

5 ^ (r * (- ?)-s (?r« (-?)> **+ ■ ■ ■ 


The convergence of all the terms in the first integral is as- 
sured when u Bt t varies as l/r* for large r , A 7 being any positive 
constant; for small r we assume u Btt varies directly with 
cross-sectional area of the body. The order of magnitude of 
J(x,8,0)/rj can then be seen to be given by 


J( z,s,0)_ A ( t 2 s 
UJ 


=o(? J lns ) 


(87) 


Since it follows from equations (85) and (56) that the magni- 
tude of <pw/UJ in the vicinity of the wing is 0 [(ts/l 2 )ln *], the 
relative error incurred by neglecting J/U m l is 0 (ts s /l A ). 

Following the method used in deriving equation (77), we 
can achieve a final simplification in the wing's perturbation 
potential. Thus, 

<Pw ~<P 2 w-\-q(%) (88) 

where 

Inn ndn (89) 

In the absence of analytical solutions for the body of revolu- 
tion problem, the evaluation of g(x) must be carried out by 
less direct methods. This will be discussed further in the 
section on applications. 


APPLICATIONS TO SEVERAL PROBLEMS INVOLVING 
SONIC FLOW 


regarded, if one prefers, as a physical interpretation of the 
result given in equation (80). 

Consider the case of compressible flow about the slender 
body of arbitrary cross section shown in figure 6. The gen- 
eral procedure is to consider first the complete three-dimen- 



Figure 6. — Schematic representation of elements of slender-body 

theory. 


In the following section, the exploitation of the results for 
sonic-flow conditions will be carried out in some detail. In 
view of the difficulty associated with transonic analysis, it 
appears likely that the equivalence relation of equation (85) 
will play an important part in the interpretation and use of 
experimental data as well as in purely theoretical predictions. 
The discussion will be concerned principally with applica- 
tions to slender wings and bodies and to the relationships be- 
tween the aerodynamic characteristics of the two configura- 
tions. It is obvious that the known basic information can be 
supplied either by theory or bv experiment and many of the 
results to be given are written with the idea that they can be 
used in this dual sense. 

RESUME OF PRINCIPAL RESULTS OF SLENDER-BODY THEORY 

It appears worthwhile, before proceeding to the examples, 
to re-examine the problem of transonic flow about slender 
bodies of arbitrary cross section from a heuristic, although 
less rigorous, point of view. This second approach may be 


sional problem and then to introduce simplifying assump- 
tions consistent with the restriction of slender plan forms. 
As in the more detailed analysis, it is again assumed that 
any shock waves are situated either entirely upstream of the 
most forward point of the body or entirely downstream of 
the most rearward point, or are normal shock waves if situ- 
ated along the length of the body. The resulting solution 
for the perturbation potential can be expressed as the sum 
of four parts. As in the preceding analysis (p% w a> (p 2 w 0 and 
(p- 2 B t are solutions of the two-dimensional Laplace equation 
as indicated. Thus ^ a corresponds to the two-dimensional 
incompressible-flow solution for translation of the cross sec- 
tion, and <p 2w t to that for the growth of the cross section. 
In addition to satisfying the prescribed boundary conditions 
at the body surface, these two terms satisfy the requirement 
that the lateral velocity components ( d<p/dy , dcp/dz) vanish 
at infinity. These terms alone do not furnish a satisfactory 
approximation, however, for cases in which S(x) is different 
from zero because <p 2w t acts like S' ( x ) In r at large r, and 
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hence d<p/dx is infinite at a large lateral distance. This error 
can be removed, however, by subtracting the term <p 2/? t cor- 
responding to the two-dimensional incompressible-flow solu- 
tion for the growth of a body of revolution having the same 
S ( x ) as the original body (thereby canceling (p 2w t at large 
r), and adding the three-dimensional solution <p B ( for flow 
about the same body of revolution. If (p n t is determined 
from linear theory, the results correspond to the familiar 
formulas of subsonic and supersonic slender-body theory. 
(This leads, in both cases, to the result given in equation (62) 
and reduces to equation (77) where, for subsonic flow, j(x) 
has the form given in equation (79a) and, for supersonic flow, 
the form given in equation (79b).) In keeping with the 
previous analysis, this function of x will be denoted by f(x) 
if it is determined from linear theory and by g(x ) if deter- 
mined from transonic theory. 

Although the linear-theory approximation is unsatisfactory 
at M oo = l, the same intuitive procedure can be extended to 
sonic flow. The desired expression follows if <p B t is deter- 
mined from the transonic differential equation. Thus, as 
given in equation (85), one has to a known order of accuracy, 

<Pw=<p 2 w ta ~hp2 w t t (90) 

and this result reduces for points near the body to 

(pw=<P2w-{-g(x) (91) 

where 

<P 2 W =<P 2 W a ~\~<p 2 w t , g(x)=lim t — ip 2h t ) (92) 

as indicated in figure 6. It is apparent that equation (90) 
has a dual basis for validity and represents either the relation 
afforded by transonic theory for M m = l, or that given by 
linear theory for other Mach numbers. The customary re- 


striction to slender wings and bodies must be observed in 
both applications. 

The power and weakness of the present intuitive reasoning 
is well illustrated by the fact that the relation given by 
equation (90) is found without recourse to the detailed inves- 
tigation of the earlier sections whereas the restriction to 
small angles of attack that enters in the simplification is 
overlooked. This deficiency stems from the fact that it is 
insufficient to assure that merely the infinite velocities be 
removed. Since the space involved is infinite, it is also nec- 
essary that certain integrals of velocity (see eq. (43)) be 
finite, and it is in connection with the attenuation of the 
velocities arising from the term <p 2w a that the deficiency 
occurs. One could, perhaps, have continued the heuristic 
reasoning but, once the principal idea has been established, 
formal analysis can be used to establish the restrictions and 
to evaluate the error terms involved. 

Some insight into the validity of the foregoing equations 
can be obtained by examining the numerical solution given 
by Yoshibara in reference 17 for sonic flow^ about a circular 
cone-cylinder at zero angle of attack. Inasmuch as it was 
not assumed either explicitly or implicitly that the per- 
turbation potential in the vicinity of the body has the form 
indicated by equation (91) (the boundary conditions were 
satisfied at the actual body surface rather than along the 
body axis), these results are particularly suited in this 
respect for the investigation of the region for which the 
simplified relation applies. On the other hand, the example 
is not ideal because the sharp corner at the shoulder violates 
the smoothness condition; it is, however, the only case for 
which a theoretical solution is available. Accordingly, 
figure 7 has been prepared so as to show the variation of 



Figure 7. — Variation of <PB,tlU m with r for several stations along the length of a cone-cylinder of semiapex 

angle 1/10 at free-stream Mach number 1. 
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<pB,JU m with r for several stations along the length of the 
cone for the case in which the semiapex angle 0 is 1/10 radian. 
Attention is called to the fact that the values for <p B t given 
in this figure are for a cone of unit length whereas the 
original values given in figure 5 of reference 17 are for a 
cone of length 10. A dotted line is also shown for each 
station representing the values obtained after subtracting 
<P 2 B JU m computed by 


<P 2 B ,t_ 1 dS j xlnr 

l ' a _ 2ir dx = T00 


(93) 


from at the same point. In order to illustrate 

further the nature of these results, figure 8 has been prepared 



Figure 8. — Variation of <pB?t — <P 2 Byt with x for several r of a cone- 

cylinder of semiapex angle 1/10 at free-stream Mach number 1. 

to show the variation of (p B , t —<P 2 Bt with x for various /*. 
The resulting values should, according to equation (92), be 
a function of x alone for small r. It can be seen from an 
examination of these figures that the difference <p B ,,—<p 2 Iit is 
indeed very nearly a function of x in most of the region for 
which results are available. Slight deviations occur in the 
immediate vicinity of the nose and at the largest distances 
from the body. The latter departures are so small, however^ 
that it is necessary to possess additional information foi 
greater distance from the body before one can determine the 
extent of the region for which the g{x) function is applicable. 

DETERMINATION OF y' (x) IN TERMS OF PRESSURE DISTRIBUTION ON A 
NONLIFTINO BODY OF REVOLUTION 

Although and hence f(x), can be calculated directly 
by means of linear theory for either distinctly subsonic or 
supersonic flow, general methods are not yet available for 
the theoretical determination of <p Btl in transonic theory. 
It is evident from its definition, however, that g(x) depends 
only on the longitudinal distribution of cross-section area 
S(x), and that its derivative can be determined from simple 
aerodynamic measurements of the flow about a slender 
nonlifting body of revolution having the same S(x) as the 
given body. From the point of view of applications, nothing 
is lost in not knowing the actual level of g(x), however, since 
knowledge of its gradient, g'(x) } is sufficient for the deter- 


mination of flow quantities such as velocity and pressure. 
Since the easiest flow quantity to measure is generally the 
pressure distribution on the surface of the body, perhaps the 
simplest way to determine g'(x ) is through a relation ex- 
pressing this quantity in terms of the pressure distribution. 
The necessary relations for the perturbation velocity poten- 
tial, <p BJ and the pressure coefficient C Vst are provided by 
equations (48) and (5), which reduce, in the vicinity of a 
slender nonlifting body of revolution, to 

^ lnr+g(x) (94) 

N 2d (dR Y_ 2 d S'\x) 

U a i>x <f>B -‘ \dx) U a dx <pB ’ ‘ 4 ttS 

(95) 

where li(x) represents the radius of the body of revolution 
and the prime denotes differentiation with respect to x. 
These relations can be combined to solve for g'(x) in terms of 
the surface pressures and the cross-section area with the 
following result: 


/(*)= 




7T 


S /2 (x ) - 1 

4tS J 


(96) 


The cone-cylinder solution of Yoshihara (ref. 17) again 
affords a means of illustrating the application of this result 
atM 00 = l. Thus, figure 9 has been prepared to illustrate 
the variations with x of the pressure coefficient 2 on the sur- 



Figure 9. — Variation of C p and g'(x) on the surface of a cone-cylinder 
of semiapex angle 1/10 at free-stream Mach number 1. 

face of a cone-cylinder having a semiapex angle of 1/10 
radian, and of g'(x)/U m computed therefrom, using equa- 
tion (96). As in the case of figures 7 and 8, the values of 
g'{x) have been converted from those given originally for a 
cone of length 10 to those for a cone of unit length. 

It is likewise evident that the function g'(x) can also be 
determined from pressure-distribution data for thin wings 
in an analogous manner, although naturally more geometric 
quantities are involved in the calculation. 


2 The curve for C PB t shown in figure 9 differs from that given originally in reference 17 
due to the correction of a sign error in the quadratic term of the expression for pressure co- 
efficient. For the cone having a semiapex angle of 1/10 radian, this change has the effect of 
diminishing the values given originally for C VB t on the cone surface by a constant amount, 
namely, 0.02. 
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RELATION BETWEEN PRESSURE DISTRIBUTIONS ON RELATED WINGS 

AND BODIES 

Wings and bodies having same longitudinal distribution 
of cross-section area. — Equation (90) displays the relation- 
ship between the perturbation velocity potential (p w for a 
thin low-aspect-ratio wing and the corresponding potential 
<p Btt for a slender nonlifting body of revolution having the 
same longitudinal distribution of cross-section area. In 
most practical applications, however, one is not so much 
interested in relations involving the velocity potential as 
those involving the pressure distributions. The following 
discussion will be concerned with the derivation of such a 
relation. Thus, consider the two objects illustrated in 
figure 10. Both have the same S(x), but the first is a non- 



Figure 10. — Views of wing and body having the same longitudinal 
distribution of cross-section area. 


lifting body of revolution and the second is a thin lifting 
wing. The relations for the potential and pressure coeffi- 
cient for the body of revolution are those given in equations 
(94) and (95) of the preceding section. The corresponding 
relations for <p w and C Vw in the vicinity of the wing are 

<Pw=<P2 w +0 0*0 (97) 


(1 p w —' 


2 c W 
U m dx 


(98) 


Since g(x) is the same for both objects, the desired relation 
between the pressure distributions on the wing and body 
of revolution can be determined by combining equations 
(95) through (98) ; thus 


Cp " ( Cp B,X^R U a ~£r + 




ST(x) 

AttS 


(99) 


It is interesting to note that this relation holds not only for 
nonlinear theory of sonic flow, but also for linearized slender- 
body theory for subsonic and supersonic flow. This follows 
directly from the fact that equation (91) and the associated 
statement are equally correct in linearized slender-body 
theory if g(x) is replaced by/(x). 

The term involving <p 2w can be considered known inasmuch 
as it can be determined directly using equation (50) or any 
of several other methods (e. g., conformal mapping, etc.) 
available from classical two-dimensional potential theory, or 
indirectly if either the linear theory or the slender-wing- 
theory solutions are known for the wing. To illustrate, let 
the subscript S denote the values indicated by the slender- 


wing-theory solution. Then, for example, if (<pw)s is 
known, <p 2[V is given by 

<P-2 W = (<Pw)s~f(x) (190) 

where /(x) is provided by equations (79). Correspondingly, 
one has 

(,0I) 

If, on the other hand, the linear-theory solution is available, 
the relation 

lim ip L =lim (p s (102) 

M oo -U Mco-U 

applies, whence 


2 

U a dx 


lim 

Moo — U 



(103) 


where the subscript L refers to values given by linear theory. 

Equation (99) enables one to calculate the pressures in 
the vicinity of any thin low-aspect-ratio wing, provided the 
pressure distribution is known on the surface of a nonlifting 
body of revolution having the same longitudinal distribution 
of cross-section area S'(x). The corresponding rule relating 
the pressures on two wings having different cross-section 
shapes but the same S(x) can be easily derived by applying 
equation (99) twice and subtracting so as to eliminate all 
terms pertinent to the body of revolution.. 

Wings and bodies having similar longitudinal distribution 
of cross-section area. — It is a simple matter to extend the 
previous results so as to include more general relations which 
apply to wings and bodies having longitudinal distributions 
of cross-section area that are merely proportional. The 
information needed to achieve this generalization is supplied 
by the transonic similarity rule for slender bodies of revolu- 
tion (ref. 5). The rule states that at M m = 1 the pressure 
distributions on two slender bodies of revolution having 
area distributions given by 



S m , ii 
S m , i 



are related according to 


(104) 


r fx\S m . u r /X\ 1 d*S u (x/l) 
** n w s m . i v »Air* d(x/iy 


In 


S mtl 

Sm, II 


th-iv* 

Tn+l / 


(105 


where S m refers to the maximum cross-section area and the 
subscripts I and II refer to properties associated with the 
two bodies. If both bodies are in air, Ti=Tn, and equation 
(105) reduces to 


a 


p B,n 



s Tr 

m,II sy 

s _ p »Ai 


1 <PSn(x/l) f S m>1 
7T d(x/iy Ln s m , u 


(106) 


If it is desired to determine the pressures C Vw for a wing 
having an area distribution given by Sn(x/l) and the pressure 
distribution is known for a body of revolution having an area 
distribution Si(x/l) proportional to Sn(x/l), one merely com- 
putes Cp BU for a body of revolution of area S u (x/l) using 
equation (106) and substitutes the result for C Vh t in equation 
(99). 
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APPLICATION TO THE CALCULATION OF PRESSURES AND FORCES ON THIN 
ELLIPTIC CONE-CYLINDERS AT ilfa>=/ 

The relations developed in the preceding section will now 
be applied to the calculation of the pressures and forces at 
A/«o = l on the conical portion of the thin cone-cylinder of 


The pressures on such a body can be determined from those 
shown graphically in figure 9 for 0=0.10 by application of 
equation (106), which reduces to the form 

^=lOO0 2 (C p J, =olo -20 2 /nlOO0 2 (112) 



elliptic cross-section illustrated in figure 11. The ordinates 
of the upper surface of the cone are given by 


z«(z,2/)=2^ (107) 


where m is the tangent of the semiapex angle, l is the length 
of the cone, and t is the maximum thickness of the cone. It 
follows that the elliptic section in the plane x=x x has major 
and minor semiaxes equal to mx x and txj2l, respectively. 
The cross-section area and surface slope are, respectively, 


cv \*tmx i 2 dz u _ mtxi 

l{mW-y 2 )' A 


(108) 


Pressure distribution on nonlifting cone-cylinders. — From 
equation (50), <p 2 w for the symmetrical nonlifting case can be 
written as 

f“ 2U„^(x, yi )ln[(y- yi y+zY*dyi (109) 


which, when evaluated 
—mx<iy<^mx), yields 




on the wing surface 

U m tmx , mx 
= - on ln~r 


(i. e., 


2 = 0 , 

(110) 


After inserting this relation into equation (99) and carrying 
out the indicated operations, one obtains 

e '*-°’.-lf( 1+,n lr) (lll) 

where C Vli refers to the pressure distribution on a circular 
cone-cylinder having such a semiapex angle 0 that it has the 
same longitudinal distribution of cross-section area as the 
elliptic cone-cylinder; thus 


As mentioned previously, following equation (99), the 
difference Cp w —C VB is the same as given by linearized 
slender-body theory for subsonic or supersonic flow. As a 
corroboration of this statement, consider the expression 
given b\ r slender-body theory for the supersonic pressure on 
the thin elliptic cone (ref. 25, p. 257) 

o.„— (1.3) 

and the corresponding expression for the supersonic pressure 
distribution on the slender circular cone (ref. 25, p. 241) 

The difference between equations (113) and (114) obviously 
reduces to the form given in equation (111). 

The application of the foregoing theory to a specific case 
will now be illustrated by determining the pressure distribu- 
tion C Pw at M m = 1 on an elliptic cone having m= % and 
£//=0.06. The first step is to calculate the pressure dis- 
tribution at M m = l on the surface of a circular cone- 
cylinder having a semiapex angle given by 

6=(^y i ==(0.015)»=0.1225 (115) 

The pressure distribution on the selected elliptic cone- 
cylinder can then be calculated through use of equation 
(111) and is 

C Pw =C Pb -0M04 (116) 

The results for both the circular and the elliptic cone- 
cylinders are shown graphically in figure 12. Note that the 
pressure distribution is independent of y in this case and 
that a single curve of C Pw versus x/l suffices to define the 
pressure on the wing. 

Drag of nonlifting cone-cylinders. — The drag D w at 

M m = 1 of thin elliptic cone-cylinders can be obtained by 
direct integration of the product of the pressure and surface 
slope and is expressible in the form 

Dw = A + p “ 2 y fj 2C »w ^ dxd y (1 17) 

V * 

where the integration is extended over the plan form and 
D e represents the contribution to the drag that results from 
a finite leading-edge radius of curvature. Since only that 
portion of <p w denoted by ^ 2 W contributes to D e , this quantity 
can be calculated in the same manner as described in refer- 
ences 25 and 28 for linear theory. Thus, the contribution 
per unit of span is, in slender theory, 


0 = 


(fr 


dD t ( *Y 

dy w 2 Tn \dx ) 


(118) 
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Figure 12. — Pressure distribution on a thin elliptic cone-cylinder at free-stream Mach number 1. 


where~r n is the radius of curvature normal to the wing 
leading edge and s is the local semispan. If the ordinate of 
the_wing, in the vicinity of the leading edge, is 

z u =h(s,y)(,s~y)M (119) 

equation (118) becomes 


dD e=I pJJJ 2 /ds\ 
dy 2 2 n {S,S) \dxJ 


( 120 ) 


For the thin elliptic cone, z u is given by equation (107) 
and s(x)=mx, hence 


and 


h 2 (s,s) 


st 2 

2m 2 / 2 


D ,= 2 mdx=-' 

Jo dy 4 


.U, 


m 2 t 2 


( 121 ) 


The second term on the right of equation (117) becomes, upon 
substitution of equation (111) for C p 


JJ 2C ^X ^ dy= - T 2 - 

w 


Lf. 


a 


dS 


” B dx 



(l +lnf)sm] 



The drag of the elliptic cone-cylinder is thus 


D w — Db 



mHHn ^ 


( 122 ) 


Note that the circular cone-cylinder and elliptic cone- 
cylinder have different values of drag at M m = 1, even though 
they have the same area distribution. As an illustration of 
the order of magnitude of the quantities involved in equa- 
tion (122), the drag at M m = l of a circular cone-cylinder 
having semiapex angle 0=0.1225, as determined by inte- 
gration of the pressure distribution shown in figure 12, is 

(£>b)*=o. m5 = 0. 00484 p J - J l 2 (123) 


whereas that of an elliptic cone-cylinder having m = 1/2 and 
</(=0.06 is 


(Pw)m. 1/2=0.00383 

l/l = 0.06 


JJ. 


u 

o'-' oc 


m 2 t 2 



Thus, although both bodies have the same area distribution, 
the drag of the elliptic cone-cylinder is less than 80 percent 
of that of the circular cone-cylinder. 

More general results for circular and elliptic cone-cylinders 
can be obtained by combining equations (123) and (124) 
with the transonic similarity rule for the drag of slender 
bodies of revolution. The latter can be derived by inte- 
gration of the corresponding relation for the pressure given 
in equation (106) and was first given by Oswatitsch and 
Berndt in reference 5. It states that the drags at M m = 1 
of two bodies of revolution having area distributions given 
by equation (104) are related assuming both bodies are in 
air, so that 7i=7n, according to 


n fS m ,n\TT, , pMJS^D-S^iO) 
Du \S m .j) l 1 1+ 2 2t 



(124) 
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For slender circular cone-cylinders S=t6 2 x 2 and equation 
(125) reduces to 

»"=(!)' [a ■ +2' 5#=-' <W** (£fj (126) 

which becomes, upon substitution of the values given by 
equation (123) for Di and 0i 


D B =—xd 4 [1.55+2 Ind 2 ] 



(127) 


The general expressions for D w and D B can be compared more 
readily if the relation 0=(raZ/2Z)^ is introduced to express 6 
in terms of in, t, and l of the thin elliptic cone having the 
same area distribution, thus 


D b = 




1.55+2 In 



(128) 


Combination of equations (128) and (122) yields the cor- 
responding result for the drag at M m = 1 of thin elliptic cone- 
cylinders 








(129) 


Before leaving the subject of similarity rules, it is of interest 
to note that equation (129) for the pressure drag of thin 
elliptic cone-cylinders is in accord with the transonic simi- 
larity rule for the pressure drag of thin wings of finite span 
(see, e. g., ref. 4). The latter is usually given in dimension- 
less form and provides that the pressure drag coefficient 
C D at M m = 1 of a family of thin nonlifting wings of affinely 
related geometry, plan-form area S p , thickness ratio r, and 
aspect ratio A satisfy the equation 

(130) 

(131) 


where 




Pad 


D u 


2 

— Sn 


and f indicates a functional dependence. If indeterminate 
forms that arise from the infinite plan-form area of a cone- 
cylinder are avoided by letting S p represent the plan-form 
area of only the conical part of the body, S p , A, and r are 
related to in, t, and l according to 


S p =ml 2 , A=4m, t=j (132) 

and equation (129) can be rewritten as follows 



It is evident from this form of the result that C+Zr is a func- 
tion of At^ alone, as required by the similarity rule. 

Pressures and forces on lifting cone-cylinders. —The 
relations summarized in figure 6 also permit the calculation 
of the pressure distribution on a thin elliptic cone-cylinder 
when inclined at a small angle of attack. To calculate this 


quantity, one must first have the expression for <p 2 x . The 
necessary expression is well known, however, since the prob- 
lem is equivalent mathematically, for the planar boundary 
conditions, to the boundary-value problem associated with 
translation of a flat plate in a two-dimensional incompres- 
sible fluid. Thus c 0 On . and on the surface of the thin 
inclined elliptic cone are 


<P 2 w.a= ±U w a{m 2 x 2 ~y 2 )^ 


U„tmx , mx . TT , , , 

<Piw= — 2/ ln~Y±L „a{m 2 x 2 — y 2 )'^ 


(134) 

(135) 


where the upper (plus) sign is to he used on the upper surface 
and the lower (minus) sign on the lower surface. After in- 
serting equation (135) into equation (99), one obtains 




2 am 2 x 

(m 2 x 2 —y 2 y^ 


(136) 


where the convention concerning upper and lower signs still 
holds. The aerodynamic loading, or the difference in pres- 
sure between the two sides of the wing, is thus 


Apw 

pJJJ/2 


<P:\-<P*wh 


4am 2 x 
{m 2 x 2 -y 2 f 


(137) 


A sketch of the load distribution is shown in figure 13. Inte- 
gration of the loading over the plan form leads to the follow- 
ing expression for lift 


£>c= P -4+ (2iram 2 l 2 ) (138) 

Although the pressure distribution at M m = 1 is not the 
same as given by linear theory, it will be recognized that the 
load distribution and lift arise solely from <&> „ and are there- 

^ ~ W ,a 



Figure 13. — -Load distribution on an inclined thin elliptic cone- 

cylinder. 
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fore the same as given by linear slender-wing theory. One 
recognizes, consequently, that the lift of any low-aspect-ratio 
wing having such a plan form that no part of the trailing edge 
extends forward of the station of maximum span is given by 

L w = P -\f- (2t ras 0 2 ) (139) 

and the drag due to lift by 

D w -D Wa _=^ L (140) 


The fraction l / 2 enters as a result of suction forces on the lead- 
ing edge. Note that the above statements also imply that all 
reciprocal and reverse flow relations of linear theory are ap- 
plicable to lifting forces at M oa =l on slender wings at small 
angles of attack. 

As in the case of drag discussed previously, equation (138) 
for the lift of thin low-aspect-ratio wings is compatible with 
the transonic similarity rule which states that the lift-curve 
slopes of a family of thin wings of finite span and affinely 
related geometry are related according to (see, e. g., ref. 4) 

t*C l =J(AtX) (141) 


where 


I*™ (.p„UJ/2)S p a 

Substitution of the geometric relations of equation 
into equation (138) for the lift, yields simply 


(142) 

(132) 


t*C l ~\At* (143) 

which is obviously in accord with the similarity rule. The 
drag due to lift given by equation (140) is in corresponding 
agreement with the appropriate transonic similarity rule. 

MOMENTUM ANALYSIS OF DRAG OF SLENDER BODIES AT A/oo = I 

The previous example of the thin elliptic cone-cylinder has 
disclosed significant differences in the drag at M m = 1 of 
elliptic and circular cone-cylinders having the same longi- 
tudinal distribution of cross-section area. Since this finding 
is contrary to the often quoted transonic area rule, it is of 
interest to study the sonic drag of a more general class of 
bodies. This will now be done using momentum methods. 

Derivation of general relation for drag. — Consider a surface 
2 which encloses a volume containing an aerodynamic body. 
The vectorial force F on the body can be determined by 
considering the pressures and flux of momentum at 2. 
In general there results 



Figure 14. — View of surface 2 used in evaluation of drag. 

surfaces normal to the x axis and situated upstream and 
downstream of the body. The remaining part of 2, denoted 
by III, is a small circular cylinder of radius R v large enough 
so that the body is entirely enclosed within the cylindrical 
surface. Since it will be assumed that the body is slender 
and smooth enough that the necessary restrictions on are 
satisfied at all stations forward of the base, but that discon- 
tinuities in geometry or velocity may occur there, the plane 
surfaces I and II will be placed at infinity upstream and at 
the base of the body, respectively. 

It is sufficient, to the order of accuracy of transonic theory, 
to approximate p and p at points near the body by 

p i MJu 
P. U m 

and 

p— 2>»=— p»-f [(1— MJ)u 2 +v 2 +w 2 ] 


Furthermore, if attention is restricted to the streamwise 
component of force, total drag of the enclosed body is given 
by 


£>=yJJ {{MJ-\)u 2 +v 2 +w 2 ] dydz-p „ J (145) 

II III 


where v T is the radial component of velocity. This expres- 
sion holds, of course, at M m = 1 and becomes 


jj ^ 2 -\-w 2 )dydz~p^ jj uVrRit 

n *iii 


zdO dx 


( 146 ) 


F=- fj (p-pjd2-ff P v[(u a +v) .ds ] (144) 

S 2 

where vector notation is used, p and p are the local static 

pressure and density, and Vis the local perturbation velocity 
vector. 

For present purposes, the surface 2 will lie taken as shown 
in figure 14. Two parts of 2 denoted as I and IT, are plane 


An alternative form for equation (146) which will prove useful 
can be obtained by replacing the surface integral over II by 
a line integral. Thus, Green’s theorem provides 

JJ jj [(!)’+(!)■] ** 

11 11 

=— j c <P JJ<pV 2 <f>dydz (147) 
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Figure 15. — View of integration contour in plane II of 2. 


where C is a curve, situated in plane II, which goes around 
the wing and also around the control surface 2, da c is an 
element of £7, and n is the unit normal drawn into the in- 
terior of £7, as illustrated in figure 15. But the relation 
<P=<p2+£/0*0 holds near the body, and therefore within C, 
hence the equation vV~0 is satisfied in II, 

J J (S+w^dy dz=-j c v da c (148) 

II 

and equation (146) for drag becomes 

D=-'~ J ' <p ^ da-p„ j j’ uv r Udedx (149) 

III 


Relation between drag of wings and bodies having the 
same area distribution. — Consider sonic flow about two 
aerodynamic shapes, one a thin wing and the other a slender 
nonlifting body of revolution, having the same longitudinal 
distribution of cross-section area. If equations (90) and 
(91) are substituted into equation (146), and if the Fourier 
expansion for <P 2 w —<p 2 B obtained from equations (58) and 
(65) is introduced into the integrals over III and the portion 
of £7 contiguous with III, one has 




Z><P2 W , ,T T dS 

d» dac+l - a (j(x)+ 


1,„ v " *•]-'- J.f 

111 in 


u B Vr B Rjdedx 


(150) 


The contour C is here divided into two parts. The inner 
portion that immediately surrounds the wing is denoted by 
C w , whereas the outer portion is denoted by Cm. It follows 
similarly that a corresponding expression can be written for 
the drag of the body of revolution. Thus 



^P2 , 

dn 


d<r c -\-l „ 


dS 

dx 


S(r,+ l.;,r to" 

'-JJ 


u B v rB Rzdd dx (151) 


where C B refers to an integration contour drawn around the 
cross section of the body. If the exterior portions of the 
control surface 2 are now selected the same for the two 
cases (i. e., surface III is the same for the wing and body), 
it follows immediately by subtraction that 


D w — D b 



dn 



dn 


- da { 


) 


(152) 


The integrals over C w can be divided into two parts after 
recalling that the analysis applies for small angles of attack 
and that <p 2yv can be written as the sum of <f> 2 w t and <p 2]Vii , 
provided the thickness does not vanish. Substitution of this 
relation into equation (152) yields 



d<p 2 

^11 .a 

dn 


d(7 c -\- 



d*V, , 

— ^ — da ( 
dn 


~L r ‘- 


(153) 


where the integrals involving the cross-product terms are 
absent since <p 2w t and its normal derivative along C w are 
even functions of 0 , whereas <p 2w a and its normal derivative 
are odd functions of z. 

Since <p 2fVa is proportional to a, it is evident that the first 
integral of equation (153) provides a contribution to the 
drag which is proportional to the square of the angle of 
attack. This quantity is exactly the vortex drag and is 
represented by the same expression at all Mach numbers. 

The difference of the two remaining integrals gives the 
difference in the drag at M m = l of a nonlifting wing and 
body having the same longitudinal distribution of cross- 
section area. Since the two integrals will not always cancel, 
the drag of the wing and body will, in general, be different. 
One can account in this way for the difference in drag of 
thin elliptic cone-cylinders and circular cone-cylinders dis- 
closed previously by integration of surface pressures. To 
show this, one must evaluate the integrals of equation (152) 
at the shoulder of the cone-cylinder ( x=l ) using the expres- 
sions for <p 2w and y? 2/i given in equations (110) and (93). In 
the integration, the contour C w extends on both sides of the 
y axis fiom — ml to -\-ml whereas the contour C B is a circle 
of radius 6l. Upon carrying out the indicated operations, 
one obtains the same result as that given previously in 
equation (122) in which the drag of the elliptic cone-cylinder 
is substantially less than that of the circular cone-cylinder. 

Special cases for which the drag of wing and body is the 
same. -Although it is important to note the difference in the 
drag of two bodies, it is perhaps even more important to 
know under what conditions the drag of the bodies is the 
same. If attention is confined to nonlifting cases so that 
<P 2 W is zero, the vortex drag vanishes, and the condition for 
the equality of the drag of a wing and body having the same 
area distribution is that the contribution of the last two 
integrals of equation (153) cancel. This condition is 
satisfied for certain large and important classes of shapes. 
One such class includes shapes that are cylindrical at the 
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base since, then, 


and 


dn 


. d(p 2 


J =°- a *" 0 

( 154 ) 

T) w =D b 

( 155 ) 


Another includes many shapes that taper to a point at the 
rear since, then, both integrals again vanish. Other classes 
for which it is more difficult to specify the geometry include 
shapes for which the integrals have equal values different 
from zero. The latter case provides some interesting situa- 
tions in which some members of a family of wings and bodies 
having the same longitudinal distribution of cross-section 
area have the same drag and others have a different drag. 
To be more specific, consider a low-aspect-ratio pointed wing 
having a straight trailing edge normal to the free-stream 
direction and smooth airfoil sections closing with a finite 
wedge angle at the rear (an example of such a wing is a 
triangular wing with biconvex profiles), and a body of revolu- 
tion having the same area distribution as the wing. Applica- 
tion of equation (153) to this pair of bodies quickly leads to 
the conclusion that the drag of the body of revolution is 
infinitely greater than that of the wing. This is apparent 
because the integral around the base of the wing is finite, 
whereas that around the base of the body is logarithmically 
infinite since 


wedge angle at the rear, each of the integrals of equation (157) 
can be written in the form 




2 U 

K a 


a (*//)_],= 


'%T(t 



( 158 ) 


Then the product At of aspect ratio and thickness ratio is the 
same for both wings, although A and r individually may be 
different. It follows, moreover, from the fact that the ratio 
z/t of the wing ordinates to the maximum thickness is the 
same function of x/f and y/s 0 for affinely related wings, that 
<p- 2yv is the same function of y/s 0 for both wings. Since s G is 
proportional to A and At is the same for both wings, it fol- 
lows that the two integrals of equation (157) have the same 
value, and both wings have the same drag. Inasmuch as it is 
only the conditions at the trailing edge that enter into the 
integrals of equation (157), similar reasoning shows that the 
two wings also have the same drag if the condition of affinely 
related geometry applies only to the cross-section shape and 
surface slopes in the x direction at the trailing edge. On the 
other hand, if the wings have merely the same longitudinal 
distribution of cross-section area, the simple relations just 
described between the various elements of equation (157) no 
longer hold, and the wings will, in general, have different 
drags. 


APPLICATION TO NONPLANAR PROBLEMS 
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and dS/dx is finite and /? is zero at x=l. The infinite drag of 
this particular body of revolution is, of course, spurious and 
is no doubt associated with the fact that the round stern is 
too blunt to treat with a theory of the slender-body type. 
On the other hand, there is no reason to believe that the 
pressure drags of the wing and body are the same. 

Since no corresponding difficulties occur at the base of the 
wing, let the drag of the above wing be compared with that 
of a second thin low-aspect-ratio wing having the same 
longitudinal area distribution. For such a comparison, 
equation (153) must be replaced by the corresponding rela- 
tion between the drag of two wings 
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It is immediately clear that the two wings have the same 
drag if they have the same geometry at the base. This con- 
dition for the equality of drag of two bodies having the same 
area distribution has been arrived at previously by Harder 
and Klunker (ref. 16) and by Berndt (ref. 29) by somewhat 
different means. As is apparent from the preceding dis- 
cussion, this condition is sufficient but not necessary. 

As a further example, consider the case where the geometry 
of the two wings is affinely related, that is, for a constant 
chord, the second wing is derived from the first by simple 
stretchings of the y and 2 dimensions. For the present class 
of wings, having straight trailing edges normal to the free- 
stream direction and airfoil sections closing with a finite 


Equation (90) expressing the relation between the per- 
turbation potential for sonic flow about a thin low-aspect- 
ratio wing and that about a slender body of revolution has 
been derived on the assumption that the boundary condi- 
tions for the wing can be specified on a planar surface. The 
development outlined in figure 6 suggests that the result can 
actually be extended to include more general classes of slender 
shapes. Accordingly, assume that equation (90) holds for 
cases involving nonplanar boundary conditions and let the 
results given in the preceding sections for the drag at M m = 1 
of thin elliptic cones be extended to include slender elliptic 
cones of any eccentricity. The analysis proceeds identically 
to that for the thin elliptic cone, the only change being that 
<P2 W must be recalculated. This is a simple problem in two- 
dimensional potential theory since <p 2w represents the poten- 
tial associated with uniform growth of an ellipse, and the 
result, when evaluated on the cone surface, is 
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Substitution of this result into equation (153) leads to the 
following relation for the drag 
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Comparison of these two expressions with the correspond- 
ing relations for thin cones given in equations (1 10) and (122) 
shows that they differ by the inclusion of an additional factor 
(l+t/2wZ) in the more general result. Although the contri- 
bution of this term is of negligible importance for cones hav- 
ing t/ml small, it is vital for nearly circular cones, and indeed 
necessary to assure the equality of D w and I) B when the ellip- 
tic cone becomes a circular cone, that is, when t/ml= 2. In 
order to illustrate this point further, figure 16 has been pre- 
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Figure 16. — Drag of elliptic cone-cylinders at free-stream Mach 

number 1 . 

pared, showing the variation with 2 ml/t, or the ratio of major 
axis to minor axis, of the drag at M m = 1 of two families of 
elliptic cones. All members of each family have the same 
longitudinal distribution of cross-section area. As indicated, 
one family is defined by mt/l= 0.02 and includes the circular 
cone-cylinder having a semiapex angle 6 of 0.10 radian, and 
the other by mt/l= 0.03 and includes the circular cone- 
cylinder having 0=0.1225 radian. The solid line indicates 
the values computed using equation (160), and the dotted 
line those computed using equation (122). In both cases, 
the drag I) B of the circular cone-cylinder is calculated from 
equation (127). As would be anticipated, the solid and 
dotted lines coincide for thin cones, but they differ consid- 
erably for circular cones (2ml/t=l). More interesting, per- 
haps, is the extent to which the drag of a family of cone- 
cylinders having the same longitudinal distribution of cross- 
section area depends on the shape of the cross section. 

The procedures applied here to the elliptic cone-cylinders 
can also be applied to many other cases, such as wing-body 
combinations, etc. For bodies having the same longitudinal 
distribution of cross-section area as a cone-cylinder one must 
merely determine the appropriate function for <p 2 w and pro- 
ceed in the same manner as for the elliptic cone-cylinders. 
For other bodies it is also necessary to have knowledge of 


either the theoretical solution or the experimental pressure 
distribution for sonic flow around a body of revolution having 
the same (or affinely related) longitudinal distribution of 
cross-section area as the given body. It should be remarked, 
however, that the extension to some of these problems in- 
volves the assumption that equation (90) applies to non- 
pi an ar cases. 

COMPARISON WITH EXPERIMENTAL RESULTS 

In the remainder of this paper, experimental data will be 
presented and a comparison made with the predictions of 
sonic slender-body theory. Although these comparisons 
may not be ideal, since experimental data for M m = 1 are 
only available for families of shapes that strain the assump- 
tions of the theory, they show remarkable agreement with 
the theory and help define the range for which the results 
may be expected to apply. 

Cone-cylinders. -The most informative class of bodies to 
investigate further with regard to comparison of theory and 
experiment is the cone-cylinder. This is because of the 
availability of not only the similarity rules, etc., but also the 
complete solution for the pressure distribution and flow 
field in the vicinity of such bodies. Experimental data are 
also available in reference 18 by Solomon for the pressure at 
several points on the surface of two rather blunt circular 
cone-cylinders at Mach numbers near unity. The tests were 
conducted on cone-cylinders having semiapex angles of 20° 
and 25° and at Mach numbers up to about 0.96. The cor- 
responding pressures at M m = l can be obtained only by 
extrapolation. The test Mach numbers are sufficiently 
high, however, that the local Mach numbers on the bod}" sur- 
face are virtually independent of the free-stream Mach 
number. These pressures are plotted in figure 17 together 
with the theoretical pressure distribution for slender cone- 
cylinders at M oa = 1. The latter were computed using 
equation (112) together with the theoretical pressure distribu- 
tion for a circular cone-cylinder having semiapex angle 6= 
0.10, see figure 9. With the exception of one point on the 
25° cone-cylinder, the theoretical and experimental values 
are in remarkable agreement, considering the bluntness of the 
cones. 




Figure 17. — Comparison of theoretical and experimental pressure distributions on two cone-cylinders at free-stream Mach number 1. 
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It would be very informative to make similar comparisons 
for cone-cylinders that are more slender, or that have non- 
circular cross sections, but the authors are unaware of any 
suitable experimental data. Studies involving bodies of 
revolution having area distributions that differ from that of 
cone-cylinders are handicapped at the present time by the 
lack of theoretical solutions for the transonic pressure dis- 
tribution, and would have to be confined to the investigation 
of such items as the range of applicability of the similarity 
rules, the existence and lateral extent of the g(x) function, 
etc. 

Wings. — Since complete solutions for sonic flow around 
three-dimensional wings have not yet been obtained, the 
following discussion must be confined to cases in which 
experimental information is known for two or more wings or 
bodies having the same or affinely related longitudinal dis- 
tributions of cross-section area. Probably the most exten- 
sive set of data of this type is that given in reference 20 for 
a large family of affinely related wings of rectangular plan 
form having NACA 63AOXX sections. Since the results 
for M m = l can be presented most concisely by using the 
variables suggested by the transonic similarity rules for 
wings of finite span (see, e. g., ref. 4), the experimental results 
for the zero-lift pressure drag and the lift-curve slope at 
= 1 are given in figure 18 by plotting C d /t h and t^C L(x as 
functions of At As shown previously by McDevitt 
(ref. 20), these data confirm the statement provided by the 
similarity rules that the results so plotted should define a 
single curve for each aerodynamic quantity. 

The curves representing the zero-lift pressure drag and 
the lift-curve slope have the same general form for high- 
aspect-ratio wings. The curves approach horizontal lines 
for the wings of larger aspect ratio, and the values for the 
lift and drag are not too different from the theoretical values 
given by Guderley and Yoshihara in references 30 and 31 
for two-dimensional sonic flow around double-wedge profiles. 
The curves approach straight lines through the origin for 
low-aspect-ratio wings and the experimental values for the 
lift-curve slope of wings having At H less than about unity 
practically coincide with the theoretical values given by 


equation (143). The corresponding theoretical values for 
the drag at Ad m = l are not known. 

Some measure of the applicability of the theoretical results 
can still be derived, however, by examining the relation be- 
tween the experimental drags of various wings. If the 
effects of the violation of the theoretical requirements that 
occurs at the round nose of the unswept leading edge of each 
of the present family of wings can be disregarded, the dis- 
cussion following equation (158) applies and all low-aspect- 
ratio wings having a given longitudinal distribution of cross- 
section area have the same drag. Inasmuch as not many 
pairs of wings of the present family actually have the same 
area distribution, a more useful statement of the result is 
that the drag is a unique function of the area distribution. 
Since the area distribution of an affinely related family of 
wings can be specified by giving, for instance, the chord l and 
the ratio S m /l 2 (or its equivalent, the product of the aspect 
ratio and the thickness ratio) of the maximum cross-section 
area to the chord squared, it follows that the drag and 
geometry of the present family of wings are related according 
to 

(pJUl 2)i 2=/ ' (~f )^' (AT) (161) 


where fi represents an unknown function of the indicated 
variables. This relation may be contrasted with that pro- 
vided by the similarity rule that states 
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where S p refers to the plan-form area. At first glance, the 
two relationships appear to bear only slight resemblance. 
It can be seen upon closer examination, however, that some 
of the apparent differences are superficial and of little or no 
significance. Thus, let equation (161) be rewritten as 
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Figure 18. — Drag and lift at free-stream Mach number 1 of a family 
of rectangular wings having NACA 63AOXX profiles. 


This appears to be the closest that the two results can be 
brought together without introducing additional restric- 
tions or approximations. Both are now concerned with the 
same quantity, C d /t **, but equation (162) states that this 
quantity is equal to some unknown function of At^, whereas 
equation (163) states that it is equal to At ^ times some 
function of At. The only way in which both results can 
be universally correct is for the functions f 2 and f 3 to be 
constants and not dependent on either At or At^. Both 
rules are not universally correct, however, since equation (163) 
is derived from transonic slender-body theory and therefore 
can be expected to apply only to wings of small aspect ratio. 
From the foregoing considerations, one can conclude that 
the drag at M m = \ of the low-aspect-ratio wings of the 
present family must depend on the geometry in such a way 
that Cd/t* is linearly proportional to At Examination 
of the drag data of figure 18 shows that the experimental 
results exhibit precisely this trend for wings of At * less 
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than about unittL An alternative interpretation of this 
result is that the drag of a number of low-aspect-ratio wings 
of the present family all having the same chord, varies as 

the square of the frontal area; that is, D Q pJUJl 2 ^ 

depends on the square of SJl 2 . Further discussion of these 
and related points appears in reference 32. 

It appears that the degree of correspondence between 
theory and experiment disclosed above for such extreme 
cases as rectangular wings of aspect ratios 3 and 4 must be 
attributed partially to the averaging influence of integration 
and that the same close correspondence may not be found 
for more detailed quantities. For example, slender-bodv 
theory indicates that the lift on low-aspect-ratio rectangular 
wings at M m = 1 is concentrated along the leading edge. 
Although pressure-distribution measurements were not in" 
eluded in the test program reported in reference 20, pitching- 
moment measurements were made from which the center-of- 
pressure position can easily be deduced. The results indi- 
cate that the center-of-pressure position at small angles of 
attack is within the first 10-percent chord at M m =l for each 
of the wings of aspect ratio 1/2 or 1, but moves progressively 
rearward for wings of larger aspect ratio. Hence the range 
of At* for which theory and experiment agree may be ex- 
pected to be less than that indicated by the integrated lift 
and drag results. On the other hand, application to wings 
of the rectangular plan form imposes severe strain on the 
slender-body assumptions, and better agreement, or a wider 
range of applicability, might be anticipated with wings of 
other plan form, such as triangular. 

Wing-body combinations. — Several comparisons between 
the experimental zero-lift drags of wing-body combinations 
and bodies of revolution having the same longitudinal dis- 
tribution of cross-section area were given by Whitcomb in 
reference 6 in connection with his discovery of the area rule- 
The bodies tested were of such geometry that the integrals 
of equation (153) are zero and the drags of wing-body com- 
binations and their equivalent bodies of revolution should 
be the same. The experimental results show excellent 
agreement in some cases, and lesser agreement in other 
cases. These results will not be discussed further here since 
the experimental data are already analyzed from the point 
of view of equality of drag in reference 6. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Apr. 2 , 1956 
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